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SUMMARY 

The  work  described  in  this  report  is  a continuation  of  an  earlier  investigation  aimed 
at  determining  the  best  turbine  configuration  for  operation  with  the  pulsating  flow  occuring 
in  a constant-volume  gas  turbine.  In  the  earlier  investigation , a radial-inflow  turbine  was 
chosen  for  analysis.  This  initial  work  is  now  extended  to  a number  of  different  two-stage, 
axial-flow  turbines. 

The  method  of  turbine  performance  prediction  for  steady-flow  operating  conditions 
was  developed  for  a single-stage,  axial-flow  turbine.  During  this  development  it  was  found 
that  existing  correlations  for  determining  turbine  loss  characteristics  at  high  negative 
angles  of  incidence  were  inadequate  and  consequently  a modified  approach  had  to  be 
developed. 

Four  basic  two-stage  turbines,  each  having  different  blade  angles,  as  well  as  several 
versions  of  one  of  these,  were  chosen  for  analysis.  The  mean  efficiencies  of  these  machines 
for  pulsating-flow  operating  conditions  were  estimated  by  using  a quasi-steady  method 
of  analysis.  The  results  obtained  for  the  selected  turbines  are  presented  and  compared  with 
the  results  for  the  radial-inflow  turbine.  A configuration  is  proposed  that  could  possibly 
give  better  mean  efficiencies  under  pulsating-flow  conditions  than  those  actually  considered. 
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ABSTRACT 

The  work  described  in  this  report  is  a continuation  of  an  earlier  investigation  aimed 
at  determining  the  best  turbine  configuration  for  operation  with  the  pulsating  flow  occuring 
in  a constant-volume  gas  turbine.  In  the  earlier  investigation,  a radial-inflow  turbine  was 
chosen  for  analysis.  This  initial  work  is  now  extended  to  a number  of  different  two-stage, 
axial-flow  turbines. 

The  method  of  turbine  performance  prediction  for  steady-flow  operating  conditions 
was  developed  for  a single-stage,  axial-flow  turbine.  During  this  development  it  was  found 
that  existing  correlations  for  determining  turbine  loss  characteristics  at  high  negative  angles 
of  incidence  were  inadequate  and  consequently  a modified  approach  had  to  be  developed. 

Four  basic  two-stage  turbines,  each  having  different  blade  angles,  as  well  as  several 
versions  of  one  of  these,  were  chosen  for  analysis.  The  mean  efficiencies  of  these  machines 
for  pulsating-flow  operating  conditions  were  estimated  by  using  a quasi-steady  method 
of  analysis.  The  results  obtained  for  the  selected  turbines  are  presented  and  compared  with 
the  results  for  the  radial-inflow  turbine.  A configuration  is  proposed  that  could  possibly 
give  better  mean  efficiencies  under  pulsating-flow  conditions  than  those  actually  considered. 
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1.  INTRODUCTION 

The  present  theoretical  work  is  a continuation  of  earlier  work  by  Erm  (1973)  aimed  at 
determining  the  best  turbine  configuration  for  operation  with  pulsating  flow  from  a constant- 
volume  combustor.  Theoretical  studies  had  indicated  that  for  small  gas  turbines,  certain  signifi- 
cant thermodynamic  advantages  might  well  be  gained  if  the  combustion  takes  place  at  constant 
volume  rather  than  at  constant  pressure.  In  the  earlier  work,  a single-stage  radial-inflow  turbine 
was  chosen  for  analysis.  Here  some  two-stage,  axial-flow  turbines  are  examined. 

Two-stage  machines  were  chosen  for  analysis  because  they  are  the  simplest  form  of  a multi- 
stage turbine  which  will  generally  be  required  in  practice.  Two  stages  would  suffice  to  indicate 
whether  variation  of  mean  flow  angles  between  early  and  later  stages  during  pulsating-flow 
operation  would  warrant  some  biasing  of  blade  angles  to  improve  the  mean  efficiency. 

The  method  of  turbine  performance  prediction  was  developed  by  examining  a single-stage 
turbine.  The  steady-flow  performance  of  this  machine  was  compared  with  the  measured  per- 
formance of  a similar  turbine;  the  comparison  indicated  that  existing  methods  for  detei mining 
turbine  loss  characteristics  at  high  negative  angles  of  incidence  were  unsatisfactory.  As  a con- 
sequence, a modified  approach  had  to  be  developed.  The  initial  analysis  was  extended  to  the 
two-stage  turbines  thereby  enabling  their  steady-flow  performance  to  be  determined.  The  mean 
efficiencies  of  these  turbines  under  pulsating-flow  operating  conditions  were  then  estimated 
using  a quasi-steady  method  of  analysis.  Further,  due  to  the  inertia  of  rotating  parts,  fluctuations 
of  rotational  speed  arising  from  the  pulsating  flow  could  be  ignored  in  this  investigation. 

The  study  has  been  concerned  with  a comparison  of  the  mean  efficiencies  of  the  turbines 
under  pulsating-flow  conditions  rather  than  with  their  work  outputs.  These  mean  efficiencies 
affect  the  performance  of  a constant-volume  gas  turbine  in  a similar  way  that  the  equivalent 
steady-flow  efficiencies  affect  the  performance  of  a more  conventional  constant-pressure  gas 
turbine. 


2.  DEVELOPMENT  OF  METHOD  OF  TURBINE  PERFORMANCE  PREDICTION  FOR 
STEADY-FLOW  OPERATING  CONDITIONS 

The  development  of  a method  of  performance  prediction  of  a single-stage  turbine  for  steady- 
flow  operating  conditions  is  described  below. 

2.1  Details  of  Turbine  Geometrical  and  Loss  Characteristics 

The  geometrical  and  loss  characteristics  of  the  turbine  are  first  considered.  The  stator 
and  rotor,  together  with  general  velocity  diagrams,  are  depicted  in  Fig.  1 in  order  to  indicate 
the  locations  of  the  various  stations  along  the  flow  path  together  with  the  sign  convention 
adopted  for  flow  and  blade  angles. 

2.1.1  Selection  of  Turbine  Blade  Characteristics  at  Reference  Diameter 

In  order  to  operate  efficiently  with  pulsating  flows,  a turbine  should  have  a high  efficiency 
over  a wide  range  of  velocity  ratios.  As  a preliminary  approach  it  was  thought  desirable  to  select 
a turbine  blading  geometry  known  to  give  a high  steady-flow,  design-point  efficiency,  together 
with  any  features  known  to  reduce  the  tendency  for  losses  to  increase  under  off-design  conditions. 

Optimisation  of  the  design-point  performance  is  assisted  by  published  parametric  studies 
by  Hawthorne  (1957),  Shaw  (1964)  and  Smith  and  Johnston  (1967).  The  first  two  of  these  studies 
are  summarised  by  Horlock  (1966).  These  studies  indicate  that  losses  in  the  blading  are  minimised 
by  choice  of  a stage  reaction  of  about  50  %,  and  stator  and  rotor  outlet  blade  angles  of  about  60°. 
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At  off-design  conditions,  an  additional  source  of  energy  loss  arises  from  the  flow  entering 
the  rotor  at  some  (positive  or  negative)  incidence  angle.  It  has  been  shown  (e.g.  Emmert  (1950)) 
that  this  loss  is  less  for  rotor  blades  with  a large  nose  radius  than  for  those  with  a small  nose 
radius,  so  that  the  former  are  to  be  preferred  for  pulsating  flow. 

Further,  the  analytical  work  of  Ainley  and  Mathieson  (1951)  showed  that  reaction-type 
rotor  blading  with  approximately  axial  inlet  blade  angles  and  relatively  high  outlet  blade  angles 
gives  less  increase  in  loss  with  incidence  angle  than  does  impulse-type  rotor  blading;  blades 
of  this  type  also  satisfy  the  above  requirements  for  optimum  design-point  performance. 

This  indicates  that  if  the  stage  loading  coefficient,  T,  is  selected  to  be  about  1 -0  then  blading 
of  this  type  should  give  a reasonable  combination  of  design-point  efficiency,  stage  loading 
and  insensitivity  to  off-design  operation.  In  order  to  select  other  blade  characteristics  the  work 
of  Ainley  and  Mathieson  (1951)  was  used  as  a guide.  The  selected  blade  characteristics  at  the 
reference  diameter  are  set  out  in  Table  I. 

A cross  section  of  a pair  of  turbine  blades  is  shown  diagramatically  in  Fig.  2 in  order  to 
indicate  turbine  blade  nomenclature.  The  selected  blade  shape  is  shown  in  Fig.  3.  The  blade 
was  chosen  to  have  a T6  profile  (Ainley  (1948))  on  a parabolic  camber  line,  which  matched  the 
selected  blade  angles,  and  tm*xlc  = 0-200.  Also  shown  in  Fig.  3 are  the  selected  design-point 
velocity  diagrams.  The  velocities,  which  are  determined  in  Section  2.1.3,  are  relatively  low  in 
order  that  the  flow  in  the  blade  passages  under  peak  pressure  conditions  remains  subsonic. 
It  is  assumed  that  both  rows  of  blading  are  designed  to  give  a 60°  outlet  flow  angle  at  the  design 
point,  and  that  this  flow  angle  is  maintained  under  all  conditions. 
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2.1.2  Analysis  and  Calculation  of  Turbine  Loss  Characteristics 

To  calculate  the  loss  characteristics  for  the  chosen  turbine  configuration,  it  is  necessary 
to  establish  the  most  appropriate  correlation  upon  which  to  base  the  losses.  Investigators  have 
proposed  axial  turbine  loss  correlations  in  terms  of  blade  geometry  and  other  factors.  Those 
of  Soderberg  (1949)  and  Ainley  and  Mathieson  (1951)  appear  to  be  the  most  authoritative: 
the  latter  takes  into  account  the  effects  of  a greater  range  of  design  and  operating  variables 
and  is  applicable  at  off-design  conditions,  whereas  the  former  correlation  is  essentially  applicable 
at  design-point  conditions  only.  Since  pulsating  flows  involve  continual  variation  of  operating 
conditions,  the  correlation  of  Ainley  and  Mathieson  was  considered  to  be  the  more  appropriate. 
However,  a revised  form  of  this  correlation,  as  proposed  by  Dunham  and  Came  (1970),  incor- 
porates several  improvements  which  result  from  the  analysis  of  more  recent  turbine  data,  and 
has  been  adopted  for  use  in  the  current  investigation. 

Verification  of  the  validity  of  the  revised  Ainley  and  Mathieson  correlation  (generally 
referred  to  as  the  revised  Ainley  correlation  in  the  remainder  of  this  report)  is  given  by  Dunham 
and  Came  (1970)  and  Dunham  and  Panton  (1973)  who  compared  predictions  with  experimental 
results.  It  should  be  realised,  however,  that  most  of  the  experimental  verification  appears  to  be 
related  to  operation  at  positive  or  low  negative  angles  of  incidence  (as  defined  in  Fig.  2),  so 
that  the  accuracy  of  this  correlation  at  large  negative  angles  of  incidence  has  not  been  subjected 
to  an  experimental  check. 

In  the  revised  Ainley  correlation,  the  blade  row  losses  are  expressed  as  a coefficient  relating 
the  decrement  in  total  pressure  to  the  dynamic  pressure  at  outlet  from  the  blade  row.  The  total 
loss  coefficient  for  a particular  blade  row,  say  123  (refer  to  Fig.  1),  is  given  by  Fris  = (Px  — P3)l 
(/»,—  p3).  This  coefficient  is  considered  to  represent  the  sum  of  component  losses  due  to  blade 
profile,  secondary  flows  and  blade  tip  clearance. 

The  profile  loss  of  a turbine  blade  is  the  loss  due  to  boundary  layer  growth  on  the  blade 
surface  and  dissipation  in  the  blade  wake,  and  is  a function  largely  of  blade  and  flow  angles. 
The  secondary  losses  are  due  to  interaction  between  the  blade  surface  flow  and  the  boundary 
layers  on  the  end  walls  and  thus  depend  mainly  on  aspect  ratio  and  blade  loading.  The  tip 
clearance  loss  is  the  increase  in  end  loss  due  to  clearance  and  is  a function  of  blade  loading 
and  the  size  and  nature  of  the  clearance. 

The  details  of  the  calculation  of  the  component  loss  coefficients  for  the  selected  blade 
geometry,  under  a wide  range  of  flow  conditions,  are  given  in  Appendix  II.  The  variation  of 
component  and  total  loss  coefficients  with  incidence  angle,  for  the  rotor,  are  shown  in  Fig.  4. 
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Superimposed  upon  this  figure  is  a curve  showing  rotor  loss  coefficients  as  determined  by  a 
modified  correlation,  in  an  attempt  to  improve  the  agreement  between  the  predicted  performance 
of  the  current  turbine  and  the  actual  performance  of  a typical  turbine  at  high  negative  incidence 
angles.  The  need  for  this  modification,  and  how  it  was  obtained,  are  set  out  in  Sections  2.3 
and  2.4  respectively. 


2.1.3  Calculation  of  Turbine  Flow  Areas 


Before  the  flow  areas  throughout  the  single-stage  turbine  could  be  calculated,  it  was  first 
necessary  to  select  appropriate  design-point  conditions  by  analysing  the  likely  flow  conditions 
within  a constant-volume  gas  turbine.  Preliminary  cycle  calculations,  as  summarised  by  Williams 
(1975),  indicated  that  a typical  engine  cycle  could  involve  pressure  pulsations  in  the  combustor 
between  3 and  9 times  atmospheric  pressure,  with  a mean  turbine  inlet  temperature  of  1200  K. 
Also,  since  an  engine  of  about  150  kW  output  is  envisaged  and  since  the  rig  test  turbine  is  likely 
to  be  operated  at  reduced  temperature,  say  600  K mean,  then  the  following  design-point  con- 
ditions were  chosen. 

ENGINE  TURBINE  RIG  TURBINE 


r,  = 1200  K 


T,  =600K 


Px  = 620  kPa  P1  = 620  kPa 

m = 0-907  kg/s  m = 1 -283  kg/s 

Afj  = 0-6  Ma  =0-6 

The  turbine  design  calculations  were  based  upon  the  expected  rig  turbine  conditions.  The 
calculations  were  performed  by  using  the  above  chosen  design-point  conditions  together  with 
the  selected  shapes  for  the  design-point  velocity  diagrams  as  given  in  Fig.  3 and  the  loss  coeffi- 
cients appropriate  to  each  blade  row  at  the  design  condition  as  given  in  Fig.  4,  and  working 
through  the  turbine  in  a step-by-step  manner. 

Flow  conditions  and  flow  areas  at  all  stations  throughout  the  turbine  are  given  in  Table  II. 


2.2  Prediction  of  Turbine  Efficiency 

The  efficiency  of  the  single-stage  turbine  was  predicted  for  a wide  range  of  operating 
conditions,  thereby  enabling  the  predicted  performance  of  this  machine  to  be  compared  with 
the  measured  performance  of  a similar,  existing  unit. 

Values  of  efficiency  can  be  determined  as  a function  of  either  inlet  total  pressure  and  outlet 
static  pressure  or  inlet  total  pressure  and  outlet  total  pressure.  The  former  specification  of 
efficiency  is  probably  the  more  appropriate  to  use  in  the  current  investigation  since  operation 
with  constant  outlet  static  pressure  is  likely  to  be  typical  of  a constant-volume  gas  turbine, 
since  efficient  outlet  diffuser  operation  under  pulsating  flow  conditions  is  very  unlikely. 

A broad  outline  of  the  thermodynamic  analysis  associated  with  the  prediction  of  turbine 
efficiencies  is  presented  in  this  section.  Details  of  the  equations  used  and  the  techniques  adopted 
are  given  in  Appendix  III. 

When  predicting  efficiencies  for  specified  operating  conditions  the  flow  conditions  at  the 
reference  diameter  are  established  by  using  conventional  one-dimensional  thermodynamic 
methods.  The  blade  row  losses  are  determined  by  using  the  loss  characteristics  given  in  Fig.  4, 
with  due  consideration  of  the  incidence  angle  at  the  rotor  inlet  for  the  particular  operating 
conditions.  A feature  of  the  analysis  is  that  it  is  quite  general  and  is  applicable  to  both  subsonic 
and  supersonic  flow  conditions  in  the  turbine. 

In  order  to  establish  the  flow  for  the  chosen  operating  conditions,  viz.  7\,  /*„  p,  and  U , 
it  is  necessary  to  adopt  an  iterative  approach.  This  process  involves  systematically  choosing 
trial  values  of  p,  and  evaluating  the  flow  variables  through  the  turbine  until  the  mass  flow  rates 
through  the  stator  and  rotor  agree  within  defined  limits.  It  is  then  possible  to  determine  the 
turbine  efficiency.  A digital  computer  has  been  programmed  to  facilitate  the  calculations. 

2.3  Comparison  of  Predicted  Efficiencies  with  Pnblislied  Data 

The  validity  of  the  method  used  to  predict  turbine  efficiencies  may  be  assessed  by  comparing 
values  of  efficiency  predicted  by  this  method  with  appropriate  published  experimental  data. 
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The  validity  of  the  revised  Ainley  correlation  for  operation  at  normal  design-point  conditions, 
i.e.  with  negligible  incidence  angle  at  the  rotor  inlet,  has  been  demonstrated  by  Dunham  and 
Came  (1970)  and  Dunham  and  Panton  (1973).  This  means  that  the  accuracy  of  prediction  of 
the  maximum  efficiency  (i.e.  the  peak  of  the  efficiency-velocity  ratio  characteristic  curve)  is  not 
now  in  question,  but  there  is  a need  in  the  present  investigation  to  check  that  the  use  of  the 
revised  Ainley  correlation  leads  to  a realistic  shape  for  the  remainder  of  the  characteristic  curve. 

Unfortunately,  there  is  very  little  published  experimental  data  covering  the  wide  range  of 
velocity  ratio  which  has  to  be  considered  in  a constant-volume  gas  turbine  application.  Panton 
(1967)  and  Dunham  and  Panton  (1973)  present  suitable  data  for  a series  of  small,  single-stage 
experimental  turbines.  None  of  the  configurations  tested  had  exactly  the  same  blade  angles  as 
those  selected  (see  Section  2.1.1)  for  the  turbine  considered  in  the  current  investigation:  their 
turbine  A had  a mean  rotor  outlet  angle  of  72°  and  turbine  C a mean  rotor  outlet  angle  of  60°, 
while  both  turbines  had  a mean  stator  outlet  angle  of  66°.  These  turbines  gave  differing  peak 
efficiencies;  that  for  turbine  C was  within  one  percentage  point  of  that  predicted  for  the  turbine 
currently  being  investigated.  However,  the  shape  of  the  efficiency-velocity  ratio  characteristic 
curve  for  both  of  these  experimental  machines  is  the  same  when  compared  on  a normalised  basis 
(i.e.  efficiency/peak  efficiency  plotted  against  velocity  ratio/velocity  ratio  for  peak  efficiency). 
Since  these  two  turbines  have  a normalised  characteristic  curve  of  the  same  shape  even  though 
geometrical  differences  exist  between  them,  it  is  reasonable  to  believe  that  the  slight  geometrical 
differences  between  the  present  turbine  and  turbines  A and  C should  not  introduce  any  signifi- 
cant variation  in  the  shape  of  the  normalised  characteristic  curve. 

The  operating  conditions  for  which  the  efficiencies  were  predicted  were  chosen  to  cover  a 
wide  range  of  velocity  ratios.  This  was  achieved  by  keeping  values  of  7",,  P,  and  p1  constant 
at  their  design-point  values  and  varying  the  value  of  V over  the  range  from  37%  to  247%  of  its 
design-point  value  in  regular  steps.  All  details  of  the  chosen  operating  conditions  together  with 
some  of  the  derived  flow  conditions  as  well  as  the  predicted  values  of  efficiency  are  shown  in 
Table  III.  The  normalised  values  of  efficiency  and  velocity  ratio  are  superimposed  upon  the 
characteristic  curve  for  turbine  A in  Fig.  5 so  that  a comparison  could  be  made  between  the 
predicted  and  experimental  results  (also  shown  on  this  figure  are  values  of  efficiency  predicted 
by  using  the  modified  correlation — see  Section  2,4). 

It  is  clear  from  the  comparison  that  the  predicted  values  of  efficiency  based  upon  the  revised 
Ainley  correlation  give  a characteristic  curve  of  the  correct  shape  at  low  velocity  ratios  (i.e. 
positive  incidence  angles  at  the  rotor  inlet)  but  have  resulted  in  optimistic  efficiencies  at  high 
velocity  ratios  (i.e.  high  negative  incidence  angles  at  the  rotor  inlet).  This  indicates  that  the  loss 
coefficients  given  by  this  correlation  for  high  negative  incidence  angles  (e  g.  —20°  to  —60°) 
are  too  low.  These  values  occur  because  the  secondary  and  tip  clearance  formulae  express  the 
losses  as  a function  of  the  overall  flow  deflection  through  the  blade  row,  and  this  tends  to  zero 
at  i = —60°  for  this  particular  blading.  In  the  practical  situation,  however,  at  this  condition, 
the  flow  probably  separates  and  deflects  on  entering  the  blade  row,  and  then  deflects  equally 
the  other  way  in  the  blade  passages:  the  extra  losses  associated  with  this  flow  condition  are 
not  taken  into  account  in  the  revised  Ainley  correlation. 

Unfortunately,  no  more  satisfactory  method  for  determining  secondary  and  tip  clearance 
loss  characteristics  is  offered  in  the  literature,  so  that  it  is  necessary  to  modify  the  total  loss 
coefficient  for  negative  angles  of  incidence. 

2.4  Modification  of  Turbine  Loss  Characteristics  at  Negative  Incidence  Angles 

If  the  alternative  loss  correlation  of  Soderberg  (1949),  as  summarised  by  Horlock  (1966),  is 
reconsidered,  then  the  data  offered  concerning  incidence  effects  is  limited  in  scope  but  suggests 
that,  for  blades  of  the  thickness  considered,  the  total  loss  is  at  its  minimum  at  about  zero  incidence 
angle,  and  may  be  greater  at  negative  incidence  angles  than  at  the  corresponding  positive  incidence 
angles.  Emmert  (1950)  showed  that  round  nosed  blades  had  slightly  greater  losses  at  negative 
incidence  angles  than  at  equivalent  positive  incidence  angles,  whereas  the  opposite  trend  was 
noticed  for  sharp  nosed  blades.  Kroon  and  Tobiasz  (1971)  analysed  experimental  data  to 
determine  the  effects  of  incidence  on  blading  efficiency,  and  showed  that  the  pressure  loss  due 
to  incidence  could  be  approximated  by  a function  of  tan*  f,  which  is  symmetrical  about  i = 0°. 

Since  the  revised  Ainley  correlation  is  believed  to  be  accurate  at  zero  and  positive  incidence 
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angles,  it  was  decided  to  retain  it  as  the  basic  correlation,  but  modify  the  total  loss  coefficient 
curve  for  the  rotor  to  the  form  Kt4«  = K.  tan*  it  + Yt46<o>  in  order  to  conform  more  closely 
with  the  loss  characteristics  noted  above.  The  modified  curve  was  fitted  to  the  revised  Ainley 
curve  at  i4  = 0°  and  ;4  = +30°,  to  give  similar  results  for  positive  angles  of  incidence.  It  is 
shown  in  Fig.  4 where  it  is  superimposed  upon  the  corresponding  revised  Ainley  curve  determined 
in  Section  2.1.2. 

In  order  to  assess  the  accuracy  of  the  modified  turbine  rotor  loss  characteristics,  it  was 
necessary  to  perform  a checking  procedure  similar  to  that  outlined  in  Section  2.3.  For  chosen 
operating  conditions,  a number  of  values  of  efficiency  were  predicted  by  using  the  efficiency- 
prediction  method  outlined  in  Section  2.2,  but  with  the  turbine  rotor  loss  characteristics  deter- 
mined by  the  modified  correlation  instead  of  the  revised  Ainley  correlation.  The  chosen  operatin'* 
conditions,  which  are  the  same  as  those  used  in  Section  2.3,  are  shown  in  Table  III  together 
with  the  predicted  values  of  efficiency.  The  normalised  values  of  efficiency  and  velocity  ratio 
are  superimposed  as  before  upon  the  characteristic  curve  for  turbine  A as  shown  in  Fig.  5 so 
that  a comparison  could  be  made  between  the  predicted  and  experimental  results.  It  is  evident 
from  the  comparison  that  the  increased  losses  at  negative  angles  of  incidence  have  depressed 
the  efficiency  at  high  velocity  ratios  further  than  was  required  to  match  the  actual  performance 
of  the  experimental  turbine. 

As  the  normalised  characteristic  curve  for  the  experimental  turbine  falls  approximately 
midway  between  the  predicted  curves  based  on  the  revised  Ainley  and  the  modified  correlations 
respectively,  it  was  decided  to  proceed  with  the  investigation  using  both  correlations  to  determine 
turbine  loss  coefficients.  In  this  way  two  estimates  of  the  effect  of  any  geometrical  change  would 
be  obtained,  the  difference  between  them  indicating  the  degree  of  uncertainty  due  to  imprecision 
in  the  determination  of  loss  coefficients,  and  the  mean  of  the  two  values  giving  some  estimate 
of  the  effect  likely  to  be  expected  in  practice.  This  was  felt  to  be  preferable  to  further  modifying 
the  total  loss  coefficient  curve  for  the  rotor  so  that  agreement  with  practical  results  was  obtained, 
which  could  lead  to  performance  assessments  being  made  without  any  idea  of  the  uncertainties 
involved. 

It  is  clear  that  a better  knowledge  of  turbine  loss  characteristics  at  high  negative  angles 
of  incidence  will  be  needed  before  accurate  optimisation  of  geometry  for  pulsating  flows  can 
be  attempted. 

3.  ANALYSIS  OF  TWO-STAGE  TURBINES  FOR  STEADY-FLOW  OPERATING 
CONDITIONS 

The  method  of  analysis  used  for  the  single-stage  turbine  was  extended  to  predict  the  per- 
formance of  the  two-stage  turbines  for  steady-flow  operation. 

3.1  Details  of  Turbine  Geometrical  and  Loss  Characteristics 

In  this  investigation  four  different  configurations  of  two-stage  turbine  were  chosen  for 
analysis. 


3.1.1  Selection  of  Turbine  Blade  Characteristics  at  Reference  Diameter 

The  four  turbines  considered  differ  only  in  blade  angles,  other  geometrical  features  being 
the  same  as  for  the  single-stage  turbine.  In  each  case,  both  stages  have  the  same  outlet  blade 
angles  and  are  designed  to  give  identical  velocity  diagrams.  The  selected  turbine  blade  charac- 
teristics at  the  reference  diameter  are  set  out  in  Table  IV. 

Build  I retained  the  blade  angles  of  the  single-stage  machine  as  a basic  configuration  while 
the  other  three  have  variations  of  blade  angle  from  this  datum  in  order  to  assess  the  effects  of 
such  changes  on  efficiency,  since  these  changes  alter  both  the  sensitivity  of  the  blading  loss 
characteristics  to  incidence  angle  variation  and  the  range  of  incidence  angles  encountered  by 
the  blading  during  a given  pulse.  These  effects  are  demonstrated  in  Fig.  6 which  depicts  the 
variation  of  the  total  loss  coefficients  with  incidence  angle  for  the  rotor  blades  of  each  turbine 
(see  Section  3.1.2),  and  also  shows  how  the  velocity  diagrams  vary  from  the  beginning  to  the 
end  of  the  standard  pulse  considered  later  in  Section  4. 
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To  simplify  the  presentation,  only  the  velocity  diagrams  appropriate  to  loss  calculations 
based  on  the  revised  Ainley  correlation  are  giv?n:  loss  calculations  based  on  the  modified  corre- 
lation lead  to  slightly  different  flow  angle  variations  during  the  cycle. 

The  range  of  incidence  angles  encountered  for  the  standard  pulse  are  indicated  on  the 
velocity  diagrams  as  well  as  on  the  corresponding  loss  curves  by  shaded  regions.  It  can  be  seen 
that  Build  2,  with  reduced  blade  outlet  angles,  experiences  a smaller  variation  in  rotor  incidence 
angle  in  the  vicinity  of  the  design  point  than  Build  1,  but  the  loss  coefficient  is  more  sensitive 
to  incidence  angle  changes.  Build  3,  with  increased  blade  outlet  angles,  has  a greater  variation 
in  rotor  incidence  angle,  but  the  loss  coefficient  is  less  sensitive  to  incidence  angle  changes. 
Build  4,  with  non-axial  rotor  inlet  angles,  is  another  configuration  with  a smaller  variation  in 
rotor  incidence  angle  but  with  a loss  coefficient  that  is  more  sensitive  to  incidence  angle  changes. 

It  was  considered  that  these  four  configurations  would  serve  to  indicate  how  blade  angles 
should  be  optimised  for  pulsating  flow. 

3.1.2  Calculation  of  Turbine  Loss  Characteristics 

The  blading  loss  characteristics  for  the  four  turbine  configurations  were  calculated  by  using 
both  the  revised  Ainley  and  the  modified  correlations. 

The  method  used  to  calculate  the  loss  characteristics  using  the  first  correlation  is  similar 
to  that  detailed  in  Appendix  II  for  the  single-stage  turbine. 

For  the  second  correlation,  relationships  of  the  form  Ft  = K tan1  / + Ktioi  were  fitted  to 
the  appropriate  revised  Ainley  curves  at  i = 0°  and  i = +30°,  as  before,  for  each  configuration 
of  rotor.  Instead  of  fitting  curves  in  this  manner  for  the  second  stators,  an  approximate  method 
was  used  in  which  the  loss  curves  for  the  rotors  were  factored  down  to  give  the  correct  loss 
coefficients  for  the  second  stators  at  »7  = 0°,  the  loss  coefficients  at  all  other  angles  of  incidence 
being  reduced  in  the  same  proportion. 

For  Build  4,  a variation  was  adopted  in  deriving  the  loss  curves  corresponding  to  the  modified 
correlation,  since  for  non-axial  blade  inlet  angles  the  limits  of  incidence  angle  are  not  symmetrical 
with  respect  to  i — 0°.  In  this  case,  different  equations  were  used  for  positive  and  negative 
incidence  angles,  i.e.  Ft  = K tan2(f  X 90/60)  -f  Ft,o>  and  Fr  = K.  tan2  (i  x 90/120)  + Ftio> 
respectively. 

The  loss  curves  derived  by  both  of  these  correlations  for  the  four  configurations  are  shown 
in  Fig.  6 for  the  first  stators  and  the  rotors  only.  The  curves  for  the  second  stators,  corresponding 
to  both  correlations,  have  been  omitted  from  this  figure  in  order  to  avoid  confusion.  They  are 
of  similar  shape  but  slightly  below  the  appropriate  rotor  curves. 

3.1.3  Calculation  of  Turbine  Flow  Areas 

For  the  two-stage  turbines,  the  design-point  values  of  Tj,  P,  and  m,  used  in  calculating 
the  turbine  flow  areas,  were  identical  to  those  used  in  Section  2.1.3  for  the  single-stage  machine. 
The  design-point  velocity  diagrams  selected  for  the  four  turbine  configurations  are  shown  in 
Fig.  6.  In  these  diagrams  two  requirements  were  satisfied;  one  was  that  the  inlet  flow  was 
at  zero  incidence  and  the  other  was  that  each  stage  of  each  two-stage  machine  gave  the  same 
work  as  the  single-stage  turbine.  Builds  1 to  3 have  a stage  loading  coefficient  of  1 0,  therefore 
the  blade  speed  for  these  three  turbines  is  the  same.  Build  4,  with  non-axial  inlet  flow,  has  a 
stage  loading  coefficient  of  2 0,  so  the  blade  speed  is  reduced  to  l/\/20  times  that  of  the  other 
turbines. 

The  design  calculations  for  each  turbine  were  performed  using  the  above  selected  design- 
point  conditions  and  velocity  diagrams  as  well  as  the  loss  coefficients  appropriate  to  each  blade 
row  at  the  design  condition.  The  method  used  is  an  extension  of  that  used  for  the  single-stage 
turbine  and  does  not  require  further  description.  Table  V summarises  the  flow  conditions  and 
flow  areas  at  all  stations  throughout  each  turbine. 

3.1.4  Redistribution  of  Flow  Angles  by  Area  Changes 

Under  pulsating-flow  conditions  it  is  found  that  the  different  blade  rows  of  each  particular 
turbine  are  subjected  to  different  mean  inlet  flow  angles.  Since  the  inlet  blade  angles  on  the 
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latter  three  blade  rows  of  each  particular  turbine  are  the  same,  the  matching  of  the  flow  to  the 
blades  throughout  the  turbines  is  not  necessarily  optimum.  Consequently  gains  in  mean  effic- 
iency may  be  possible  if  the  matching  throughout  the  turbines  is  improved. 

In  an  attempt  to  obtain  improved  mean  efficiencies,  the  flow  areas  through  the  turbine 
were  changed,  thereby  redistributing  the  flow  directions  and  thus  altering  the  matching.  This 
procedure  was  relatively  simple  and  involved  only  minor  changes  to  the  efficiency-prediction 
programme  (referred  to  in  Section  3.2).  In  the  current  investigation  only  the  flow  areas  of  Build  1 
have  been  changed.  Altogether  seven  different  versions  of  this  build  were  investigated.  They 
were  obtained  by  decreasing  (or  increasing)  the  flow  areas  at  the  turbine  inlets  and  corres- 
pondingly increasing  (or  decreasing)  the  flow  areas  at  the  turbine  outlets  and  changing  the  flow 
areas  at  the  intermediate  locations  by  proportional  amounts.  The  details  of  the  different  versions 
are  shown  below. 


Change  of  flow 

Change  of  flow 

areas  at  turbine 

areas  at  turbine 

inlet  (station  1) 
o/ 

outlet  (station  12) 

Version  1 

/o 

0 

/o 

0 

Version  2 

-10 

+10 

Version  3 

-20 

+20 

Version  4 

- 5 

+ 5 

Version  5 

-15 

+ 15 

Version  6 

+ 10 

-10 

Version  7 

+20 

-20 

It  is  to  be  noted  that  due  to  the  bulkiness  of  the  complete  set  of  results  associated  with  all  of 
the  versions,  only  a representative  sample  is  presented  here. 

3.2  Prediction  of  Turbine  Efficiency 

The  method  used  to  predict  the  efficiencies  of  the  two-stage  turbines  for  specified  operating 
conditions,  T„  P„  pl3  and  U,  is  an  extension  of  that  outlined  in  Section  2.2  and  Appendix  III 
for  the  single-stage  turbine. 

The  iterative  process  used  to  establish  the  flow  conditions  throughout  these  turbines  is 
slightly  more  complex  than  that  used  for  the  single-stage  machine  since  in  this  case  three  pressures 
have  to  be  determined  by  iteration.  For  an  arbitrarily  selected  value  of  p7  a series  of  values 
of  p3  are  systematically  chosen  until  the  mass  flow  rates  through  the  first  stator  and  first  rotor 
agree  within  defined  limits  and  also,  for  the  same  value  of  p7,  a series  of  values  of  p,  are  system- 
atically chosen  until  the  mass  flow  rates  through  the  second  stator  and  second  rotor  agree 
within  similar  limits.  The  mass  flow  rates  through  the  first  and  second  stages  are  then  compared 
and  if  there  is  a large  discrepancy  then  the  value  of  p,  is  modified  and  the  entire  procedure 
repeated  until  the  mass  flow  rates  through  both  stages  agree  within  the  defined  limits.  When 
this  occurs,  all  other  flow  variables  form  a balanced  solution  for  the  chosen  operating  conditions 
and  it  is  then  possible  to  determine  the  turbine  efficiency. 

3.3  Prediction  of  Steady-Flow  Efficiencies  for  Operating  Conditions  Encountered 
Momentarily  by  Turbines  During  Pulsating  Flow 

When  the  turbines  are  operating  under  pulsating-flow  conditions  their  efficiencies  are 
continually  changing  because  of  continual  changes  in  the  operating  conditions.  The  efficiency- 
prediction  method  outlined  in  the  previous  section  was  used  to  predict  the  efficiencies  of  the 
two-stage  turbines  for  a number  of  discrete  sets  of  operating  conditions  encountered  by  the 
turbines  during  pulsating-flow  operation. 

A representative  sample  of  the  chosen  operating  conditions  for  which  the  efficiencies  were 
predicted  is  given  in  Table  VI  for  Build  1.  The  inlet  temperatures  and  pressures  are  typical 
combinations  obtained  from  a preliminary  analysis  of  a turbine  operating  under  a selected 
type  of  pulsating  flow  (see  Section  4.2  and  Fig.  9).  The  blade  speeds  correspond  to  the  design 
speed  and  the  outlet  pressures  correspond  to  atmospheric  pressure.  The  predicted  values  of 
efficiency  together  with  corresponding  values  of  velocity  ratio  are  also  given  in  Table  VI.  The 
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efficiencies  and  velocity  ratios  have  been  computed  as  a function  of  inlet  total  pressure  and 
outlet  total  pressure  and  also  inlet  total  pressure  and  outlet  static  pressure  for  the  cases  when 
the  turbine  loss  characteristics  are  determined  by  both  the  revised  Ainley  and  the  modified 
correlations.  The  data  contained  in  Table  VI  are  shown  plotted  in  Fig.  7 together  with  similar 
plots  for  the  other  turbines. 

In  order  to  establish  simplified  mathematical  relationships  between  r,  and  v for  use  in  the 
pulsating-flow  analysis,  and  thus  avoid  frequent  resort  to  the  detailed  turbine  flow  analysis 
described  in  Section  3.2,  the  curves  shown  in  Fig.  7,  which  correspond  to  typical  constant- 
volume  operating  conditions,  were  fitted  by  the  method  of  least  squares. 

To  be  able  to  select  an  optimum  blade  speed  for  each  turbine,  the  above  procedure  was 
repeated  for  two  additional  speeds,  all  other  operating  conditions  remaining  unchanged.  The 
additional  speeds  were  25%  above  and  25%  below  the  design  speeds. 

The  procedure  described  above  for  Builds  1 to  4 was  extended  to  the  different  versions  of 
Build  1 . The  operating  conditions,  efficiencies  and  velocity  ratios,  as  well  as  the  efficiency-velocity 
ratio  curves,  associated  with  the  different  versions  of  Build  1 are  not  documented,  however, 
since  they  are  basically  the  same  as  those  already  presented. 

4.  ESTIMATION  OF  TURBINE  MEAN  EFFICIENCY  UNDER  PULSATING-FLOW 
CONDITIONS 

In  the  preceding  work  the  two-stage  turbines  have  only  been  analysed  for  steady-flow 
operating  conditions.  The  analysis  is  now  extended  to  the  more  complex  case  of  pulsating  flow. 
Estimates  are  made  of  the  mean  efficiencies  of  the  different  two-stage  turbines  when  they  are 
operating  under  pulsating-flow  conditions  similar  to  those  that  occur  in  a constant-volume  gas 
turbine. 

The  unsteady  flow  used  in  the  current  investigation  is  defined  primarily  in  terms  of  a form 
of  cyclic  pressure  fluctuation  in  the  combustor  (i.e.  at  the  turbine  inlet).  The  form  of  the  idealised 
pressure  pulse  chosen  for  analysis  is  depicted  in  Fig.  8 for  a complete  cycle.  As  can  be  seen  from 
this  figure,  the  pulsed  cycle  is  conveniently  divided  into  three  distinct  parts,  namely  pressure 
rise  (combustion),  blowdown  and  scavenging.  Instantaneous  pressure  rise  due  to  combustion  is 
depicted  in  the  figure  although  for  more  accurate  simulation  extra  time  for  combustion  could 
be  added  to  the  cycle  time.  During  the  blowdown  phase  there  is  no  cold  scavenge  air  flowing 
into  the  combustor  from  the  compressor  and  the  pressure  continually  decreases  with  time,  as 
shown,  as  the  hot  gas  is  discharged  through  the  turbine.  During  scavenging  the  cold  incoming 
air  is  assumed  to  intermix  completely  with  the  hot,  undischarged  combustor  contents  and  the 
pressure  in  the  combustor  remains  essentially  constant. 

In  order  to  define  the  pulsating  flow  more  precisely  it  is  necessary  to  specify  certain  basic 
cycle  features.  In  an  earlier  study  on  a radial-inflow  turbine  by  Erm  (1973),  a pulse  corresponding 
to  a compressor  pressure  ratio  of  2-85:1  and  a combustion  pressure  ratio  of  3-00:1,  making 
a peak  pressure  ratio  of  8-55  :l,  was  used.  Since  this  pressure  ratio  would  result  in  choking  in 
a two-stage  turbine  (probably  three  stages  would  be  needed  in  practice  for  this  pressure  ratio), 
a pulse  corresponding  to  a compressor  pressure  ratio  of  1-71  :1  and  a combustion  pressure 
ratio  of  2-34  :1,  giving  a peak  pressure  ratio  of  4-00  :1,  was  used.  These  compressor  and  com- 
bustion pressure  ratios  were  chosen  since  they  led  to  approximately  the  same  ratio  of  velocity 
ratio  for  the  two-stage  turbines,  i.e.  the  final  value  is  2-34  times  the  initial  value  (expressed  as 
a function  of  inlet  total  pressure  and  outlet  static  pressure),  as  for  the  radial-inflow  turbine 
mentioned  above,  thus  maintaining  a comparable  severity  of  pulsation.  Variables  associated 
with  the  pulse  depicted  in  Fig.  8 as  well  as  the  turbines'  speeds  are  listed  below. 
ptv  = 173-3  kPa 

Pow  = 405  -6  kPa 

p„  = 101  -4  kPa 

U = 184-8,  246-5,  308-1  m/s  (Builds  1 to  3)  (Build  1,  Versions  1 to  7) 

U = 130-7,  174-3,  217-8  m/s  (Build  4) 

fV av  = 173-3  kPa 

Tscav  = 359-1  K 

The  reasoning  behind  the  selection  of  values  for  pov  and  />ow,  the  pressures  in  the  combustor 
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before  and  after  combustion  respectively,  has  already  been  given.  The  values  of  U correspond 
to  the  design  speeds  and  speeds  25%  above  and  below  the  design  speeds.  The  values  of  Pscav 
and  Tscav  correspond  to  the  pressure  and  temperature  respectively  of  the  scavenging  air  after 
the  air  has  been  compressed  from  atmospheric  conditions  through  a pressure  ratio  of  1 -80 :1 
and  has  undergone  a charging  pressure  loss  of  5%.  The  combustor  volume,  v0,  was  chosen  to 
be  0-0566  m*  (2 -Oft.*),  the  same  as  that  used  for  the  radial-inflow  turbine.  The  values  of  total 
temperature  and  pressure  at  the  turbine  inlet  are  taken  to  be  equal  to  their  respective  static 
values  within  the  combustor. 

As  was  done  for  the  radial-inflow  turbine,  the  pulsating  flow  in  the  current  investigation 
was  analysed  by  using  a quasi-steady  approach.  In  order  to  do  this  it  was  necessary  to  sub- 
divide the  pulsed  cycle  into  a number  of  intervals,  each  corresponding  to  a time  increment,  Ar, 
in  the  manner  shown  in  Fig.  8.  Uniform  intervals  need  not  necessarily  be  used,  but  for  simplicity 
they  were  used  here.  It  was  assumed  that  the  combustor  pressure  and  temperature  remained 
constant  for  the  duration  of  these  short  time  intervals.  The  justification  for  analysing  the  pulsating 
flow  in  this  manner  is  discussed  in  Section  5.2.2. 

The  pulsating-flow  analysis  is  simplified  by  assuming  that  all  flow  is  one-dimensional  and 
that  the  specific  heats  of  the  gas  do  not  vary  with  temperature.  In  order  to  simplify  the  presen- 
tation, only  a broad  outline  of  the  analysis  is  given  in  this  section.  Further  details  are  given  in 
Appendix  IV.  Both  the  blowdown  and  scavenging  phases  of  the  pulsating  flow  are  analysed 
thermodynamically  and  then  the  method  used  to  estimate  the  mean  efficiency  of  a turbine  when 
it  is  operating  under  pulsating-flow  conditions  is  described. 

4.1  Blowdown  Phase 

Before  the  blowdown  calculations  can  be  commenced,  it  is  first  necessary  to  select  a trial 
value  of  /ov,  the  combustor  temperature  before  combustion.  It  will  generally  be  necessary  to 
modify  this  value  of  fov  at  a later  stage  if  it  does  not  agree  with  the  computed  combustor  tempera- 
ture at  the  end  of  the  scavenging  phase. 

The  mass  of  gas  in  the  combustor  at  the  start  of  blowdown  can  be  determined  from  the 
initial  pressure,  temperature  and  volume.  A pseudo  mass  flow  rate  can  be  determined  by  using 
these  initial  conditions  (see  Appendix  IV  for  details)  and  from  this  mass  flow  rate  it  is  possible 
to  determine  the  mass  outflow  during  the  first  time  interval.  The  actual  increment  of  turbine 
work  in  the  interval  can  be  determined  by  first  computing  the  isentropic  increment  of  work 
from  a knowledge  of  the  increment  of  mass  flow,  the  gas  temperature  and  pressures,  and  then 
multiplying  the  result  by  the  turbine  efficiency.  The  efficiency  can  be  determined  by  computing 
the  velocity  ratio,  which  depends  upon  the  initial  temperature  and  pressures  as  well  as  the  speed, 
and  then  using  the  appropriate  efficiency-velocity  ratio  relationship  shown  in  Fig.  7. 

At  the  end  of  the  time  interval,  the  residual  mass  of  gas  in  the  combustor  can  be  determined 
and  the  new  pressure  and  temperature  calculated  assuming  that  the  combustor  contents  have 
expanded  isentropically  to  the  new  condition.  These  new  conditions  become  the  initial  conditions 
for  the  following  interval.  The  procedure  can  be  repeated  for  successive  time  intervals  until 
the  combustor  pressure  falls  to,  or  below,  the  scavenging  pressure,  Pscav,  at  which  point  the 
scavenging  process  is  commenced. 

Throughout  the  blowdown  process  the  isentropic  work  increments  and  the  actual  turbine 
work  increments  are  summated,  with  the  view  of  obtaining  totals  for  the  whole  cycle  from 
which  the  mean  efficiency  can  be  calculated. 

4.2  Scavenging  Phase 

The  scavenging  calculations  are  very  similar  to  those  outlined  above  for  the  blowdown 
phase.  In  this  case,  however,  an  allowance  has  to  be  made  for  the  effects  of  the  incoming 
scavenging  air. 

For  the  first  time  interval  of  the  scavenging  phase  the  turbine  entry  conditions  are  equal 
to  those  at  the  end  of  blowdown  and  the  mass  outflow  and  the  isentropic  and  actual  increments 
of  work  can  be  determined  by  using  the  method  outlined  above  for  the  blowdown  phase.  During 
the  interval,  the  volume  outflow  is  replaced  by  an  equal  volume  of  entering  scavenging  air. 
At  the  end  of  the  time  interval,  the  mass  of  gas  in  the  combustor  can  be  determined  and  the 
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temperature  of  the  mixture  calculated  on  the  assumption  that  the  cold  incoming  air  intermixes 
completely  with  the  hot  undischarged  combustor  contents.  This  new  temperature  becomes  the 
initial  temperature  for  the  following  interval. 

The  procedure  can  be  repeated  for  successive  time  intervals  until  scavenging  is  terminated 
at  the  point  where  a specified  mass  of  cold  scavenging  air  has  entered  the  combustor.  This  mass 
was  chosen  to  be  equal  to  the  mass  contained  within  the  combustor  at  scavenging  pressure 
and  temperature  (as  for  the  radial  turbine). 

Throughout  the  scavenging  process  the  isentropic  work  increments  and  the  actual  work 
increments  are  summated  and  added  to  those  for  the  blowdown  phase  to  obtain  totals 
for  the  cycle. 

At  the  end  of  the  scavenging  phase,  a comparison  is  made  between  the  final  value  of  /„ 
the  combustor  temperature,  and  the  trial  value  of  rov  selected  at  the  start  of  the  blowdown 
analysis.  If  there  is  a large  discrepancy  between  these  two  values  of  temperature,  then  the  value 
of  /#v  is  modified  and  the  cycle  calculations  repeated  until  the  agreement  is  acceptable. 

Table  VII  lists  a number  of  flow  characteristics  at  various  stages  throughout  a pulsating-flow 
cycle.  The  data  presented  in  Table  VII  are  only  a representative  sample  of  the  results  obtained, 
and  correspond  to  Build  1 when  it  is  operating  at  the  design  speed  and  when  the  turbine  loss 
characteristics  are  determined  by  the  modified  correlation.  The  combustor  pressure,  temperature 
and  mass  listed  in  Table  VII  are  shown  plotted  in  Fig.  9.  The  twenty  three  pressure-temperature 
combinations  which  were  used  when  predicting  turbine  efficiencies  (see  Section  3.3)  are  clearly 
indicated  on  this  figure.  A number  of  corresponding  turbine  internal  flow  characteristics,  most 
of  which  are  not  listed  in  Table  VII,  are  shown  plotted  in  Fig.  10.  Detailed  discussion  on  this 
figure  is  given  in  Section  S. 


4.3  Estimation  of  a Mean  Turbine  Efficiency 

Now  that  both  the  blowdown  and  scavenging  phases  of  the  pulsed  cycle  have  been  analysed 
thermodynamically  it  is  possible  to  estimate  the  mean  total-static  efficiencies  of  the  turbines 
when  they  are  operating  under  pulsating  flow.  The  mean  efficiency  is  simply  the  ratio  of  the  sum 
of  the  actual  work  increments  to  the  sum  of  the  isentropic  work  increments  for  the  whole  cycle. 

The  mean  total-static  efficiencies  of  Builds  1 to  4 when  they  are  operating  at  75%,  100% 
and  125%  of  ther  design  speeds  and  when  their  loss  characteristics  are  determined  by  both 
the  revised  Ainley  and  the  modified  correlations,  are  shown  in  Fig.  1 1 . The  maximum  mean 
total-static  efficiency  that  could  reasonably  be  expected  in  practice  for  the  pulsating  flow 
considered  is  approximately  75-6%  and  would  be  obtained  by  Build  3 when  operating  at  about 
its  design  speed. 

The  mean  total-total  and  the  mean  total-static  efficiencies  of  Build  1,  Versions  1,  3 and  6, 
when  they  are  operating  at  the  design  speed  of  Build  1 and  when  their  loss  characteristics  are 
determined  by  the  modified  correlation,  are  shown  in  Fig.  12. 

Discussion  on  the  above  results  is  given  in  Section  5. 


5.  DISCUSSION 

A method  of  turbine  performance  prediction  for  steady-flow  operating  conditions  has  been 
developed  initially  for  a single-stage  turbine  and  extended  to  a number  of  different  two-stage 
turbines  thereby  enabling  their  steady-flow  performance  characteristics  to  be  determined.  Once 
the  steady-flow  behaviour  of  these  turbines  is  known,  then  their  mean  efficiencies  when  operating 
under  a selected  type  of  pulsating  flow  can  be  estimated. 

In  this  work  a number  of  factors  arose  which  require  examination  in  order  to  assess  their 
effects  upon  the  validity  of  the  conclusions.  The  discussion  of  these  factors  is  grouped  into 
sections  which  deal  with  the  steady-flow  ard  pulsating-flow  operation  of  the  turbines  respectively 
and  also  a section  that  compares  the  estimated  mean  efficiencies  of  the  axial-flow  turbines  with 
those  of  a radial-inflow  turbine  (Erm  (1973)). 
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5.1  Steady-Flow  Operating  Conditions 

Various  aspects  of  the  analysis  associated  with  the  steady-flow  operation  of  the  turbines 
are  now  considered. 

5.1.1  Turbine  Loss  Characteristics 

When  the  theoretically  predicted  values  of  efficiency  for  the  single-stage  turbine  were 
compared  with  corresponding  data  of  Panton  (1967)  significant  differences  between  the  two 
sets  of  results  for  operation  at  high  negative  angles  of  incidence  suggested  that  the  revised  Ainiey 
correlation  is  inadequate  to  determine  the  turbine  loss  characteristics  in  this  region  of  operation. 
One  possible  explanation  for  these  differences  in  efficiency  is  that  the  turbines  used  in  the  com- 
parison are  not  identical.  However,  this  explanation  can  be  discounted  since  it  was  shown  in 
Section  2.3  that  two  turbines,  somewhat  similar  to  the  single-stage  turbine  considered  in  this 
investigation,  had  a normalised  efficiency-velocity  ratio  characteristic  curve  of  the  same  shape 
even  though  definite  geometrical  differences  existed  between  them. 

In  an  attempt  to  obtain  better  agreement  between  theoretical  and  experimental  results,  the 
rotor  loss  characteristics  of  the  single-stage  turbine  were  modified  to  theform  Vt  = Ktan*/  + yT«». 
Justification  for  modifying  the  turbine  loss  characteristics  in  this  manner  is  given  in  Section  2.4. 
A comparison  of  predicted  and  actual  performance  showed  that  the  modified  correlation  was 
also  inadequate  to  determine  the  turbine  loss  characteristics  at  high  negative  angles  of  incidence; 
in  fact  the  use  of  the  modified  correlation  overcorrected  the  efficiency  disparity  and  led  to 
efficiencies  as  far  removed  from  the  experimental  normalised  characteristic  curve  as  initially. 

It  is  clear  from  the  above  that  the  exact  form  of  the  loss  characteristics  of  the  single-stage 
turbine  considered  in  this  investigation  at  high  negative  angles  of  incidence  is  not  known; 
however,  it  could  be  expected  that  the  losses  are  about  midway  between  those  determined  by 
using  both  correlations.  Consequently  it  was  decided  to  proceed  with  the  planned  investigation 
using  both  correlations  since  by  doing  this  at  least  the  limits  would  be  obtained  of  the  true 
turbine  behaviour  to  be  expected  in  practice.  Some  justification  for  doing  this  is  given  in 
Section  2.4.  It  is  to  be  emphasised  that  very  little  experimental  investigation  of  turbine  operation 
at  high  negative  angles  of  incidence  has  been  reported,  and  therefore  the  exact  nature  of  turbine 
loss  characteristics  in  this  region  of  operation  is  uncertain.  More  data  are  needed  before  reliable 
optimisation  of  turbine  geometry  for  pulsating-flow  operation  can  be  attempted. 

5.1.2  Prediction  of  Turbine  Efficiency 

The  use  of  a simplified  analysis  in  the  current  investigation  raises  the  question  whether 
or  not  the  results  obtained  are  reliable.  A simplified  analysis,  similar  to  that  used  in  this  investig- 
ation, was  used  by  Dunham  and  Came  (1970),  who  compared  predicted  values  of  efficiency  of 
twenty  five  turbines  with  corresponding  experimentally  measured  values.  The  predicted  values 
of  efficiency  were  within  2%  of  the  corresponding  experimentally  measured  values  except  in 
a few  cases.  The  good  agreement  between  predicted  and  measured  efficiencies  obtained  by  these 
workers  by  using  a simplified  thermodynamic  analysis  gives  confidence  in  the  accuracy  of  the 
predicted  values  of  efficiency  obtained  in  the  current  investigation. 

5.2  Pulsating- Flow  Operating  Conditions 

Various  aspects  of  the  analysis  associated  with  the  pulsating-flow  operation  of  the  turbines 
are  now  considered. 

5.2.1  Idealised  Flow  into  Turbine 

The  flow  into  the  turbine  is  idealised  and  is  based  upon  a number  of  approximations  which 
have  been  discussed. 


The  use  of  this  idealised  pulsating  flow  provides  a simple  means  of  making  a comparison 
between  the  unsteady  and  steady  performances  of  a particular  turbine.  Similarly,  it  is  possible 
to  compare  the  unsteady  performance  of  different  turbines. 


5.2.2  Quasi-Steady  Method  of  Analysis 

The  pulsating  flow  used  in  the  current  investigation  was  analysed  by  using  a quasi-steady 
method.  This  means  that  the  unsteady  flow  conditions  existing  throughout  the  turbines  were 
approximated  by  a series  of  steady-flow  processes.  Although  not  as  desirable  as  a more  rigorous 
treatment,  this  was  a tractable  way  of  dealing  with  the  physics  of  pulsating  flow.  The  question 
does  arise,  however,  of  the  accuracy  of  this  approach. 

A justification  for  use  of  the  quasi-steady  method  of  analysis  is  that  the  rate  of  change  of 
the  turbine  inlet  conditions  is  slow  enough  that  the  inlet  conditions  do  not  change  significantly 
during  the  time  required  for  a particle  of  gas  to  pass  through  the  turbines.  For  the  two-stage 
turbines  operating  at  optimum  velocity  ratio,  with  typical  temperatures  and  pressures,  this  time 
is  of  the  order  of  1 ms. 

Furthermore,  as  indicated  by  Horlock  (1968),  the  value  of  the  frequency  parameter,  given 
by  u = (frequency  x blade  chord)/(2  x mean  velocity),  can  be  used  to  assess  the  accuracy 
of  a quasi-steady  method  of  analysis.  For  the  present  investigation  the  value  of  the  frequency 
parameter  for  the  four  turbines  is  in  the  region  of  0-01.  Horlock  (1968)  indicates  that  quasi- 
steady analyses  appear  to  be  of  value  at  these  values  of  frequency  parameter. 

It  is  also  essential  that  the  time  intervals  used  in  the  calculations  are  kept  small  enough 
that  the  turbine  inlet  conditions  do  not  change  excessively  from  one  time  interval  to  the  next. 
In  the  current  investigation,  the  time  intervals  were  chosen  to  be  0-001  s.  Table  VII  lists  selected 
data  at  0-010  s intervals. 

It  is  apparent  from  Table  VII  and  Fig.  9 that  the  above  two  requirements  are  satisfied.  Any 
inaccuracies  resulting  from  the  use  of  this  method  will  primarily  be  due  to  the  regions  of  more 
rapid  temperature  and  pressure  change.  The  possible  inaccuracies  will  be  minimised,  however, 
due  to  the  fact  that  these  regions  are  only  a small  portion  of  the  total  cycle  period. 

5.2.3  Turbine  Configuration  Giving  Best  Mean  Efficiency 

The  mean  total-static  efficiencies  of  the  four  basic  turbines  for  three  speed  ratios  and  for 
losses  determined  by  two  correlations  are  shown  in  Fig.  1 1 . Curves  have  been  drawn  through 
the  mean  efficiencies  of  each  turbine  for  each  loss  correlation.  Due  to  the  fact  that  only  three 
points  were  available  for  each  curve,  then  the  curves  have  been  simplified  to  arcs  of  circles.  The 
two  individual  mean  efficiency  curves  for  each  particular  turbine  have  been  averaged  as  shown 
in  order  to  obtain  additional  curves  that  correspond  to  an  average  mean  efficiency.  If  it  is 
assumed  that  the  averaged  mean  efficiencies  are  the  mean  efficiencies  that  could  reasonably  be 
expected  in  practice,  then  some  tentative  conclusions  can  be  drawn.  Before  this  is  done,  however, 
it  is  convenient  to  compile  the  following  list  (approximate  values  only). 


Maximum 

Speed  corresponding  to 

averaged 

maximum  averaged  mean 

mean 

efficiency  expressed  as  a 

efficiency 

0/ 

percentage  of  design  speed 

*»/ 

Build  1 

/o 

74-7 

/ O 

100-8 

Build  2 

70-1 

97-5 

Build  3 

75-6 

100-0 

Build  4 

70-7 

115-8 

Clearly  the  performance  of  Build  3 is  superior  to  that  of  the  other  turbines  for  the  entire 
speed  range  considered  and  this  turbine  would  be  the  best  to  use  for  the  pulsating  flow  considered 
in  this  investigation.  The  performance  of  Build  1 is  always  better  than  that  of  Builds  2 and  4. 
Although  the  best  performance  of  Build  4 is  better  than  the  best  performance  of  Build  2,  the 
choice  between  these  two  turbines  depends  upon  the  speed  range  considered.  At  low  speeds 
Build  2 is  preferable  to  Build  4 while  at  high  speeds  the  converse  is  true. 

An  analysis  of  Fig.  7,  which  depicts  the  steady-flow  efficiency-velocity  ratio  curves  for  the 
four  turbines,  indicates  that  Builds  1 and  2,  with  outlet  blade  angles  of  60°  and  50°  respectively, 
produce  the  best  total-total  efficiencies  and  Builds  1 and  3,  with  outlet  blade  angles  of  60°  and 
70°  respectively,  produce  the  best  total-static  efficiencies.  This  suggests  that  a turbine  with  outlet 
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blade  angles  in  the  range  50°  to  60°  in  the  first  stage  and  outlet  blade  angles  in  the  range  60° 
to  70°  in  the  second  stage  would  give  the  best  total-static  efficiencies  under  steady-flow  conditions 
and  better  mean  total-static  efficiencies  under  pulsating-flow  conditions  than  the  four  turbines 
investigated. 

S.2.4  Redistribution  of  Flow  Angles  by  Area  Changes 

An  analysis  of  Fig.  10,  which  depicts  the  manner  in  which  the  internal  flow  characteristics 
vary  throughout  Build  1 when  it  is  operating  at  its  design  speed,  indicates  that  the  incidence 
angles  on  the  second  rotor  vary  more  than  those  on  the  first,  and  are  biased  toward  the  negative 
side  over  most  of  the  cycle.  Also,  the  second  stage  does  the  most  work  early  in  the  cycle  but 
does  the  least  work  later.  Consequently,  any  change  of  the  second  rotor  inlet  angle  in  the  negative 
direction  to  give  a better  average  alignment  to  the  flow  could  be  detrimental  in  the  early  part 
of  the  cycle,  when  the  turbine  is  highly  loaded,  but  to  advantage  over  the  rest  of  the  cycle. 

An  attempt  to  obtain  improved  mean  efficiencies  was  made  by  changing  the  flow  areas 
throughout  the  turbine  thereby  redistributing  the  flow  directions.  The  details  of  the  different 
versions  of  Build  1 analysed  in  this  investigation  are  given  in  Section  3.1.4.  This  method  of 
altering  the  matching  between  the  gas  flow  angles  and  the  inlet  blade  angles  is  not  as  desirable 
as  modifying  the  inlet  blade  angles,  but  the  method  was  relatively  simple  and  involved  only 
minor  changes  to  the  efficiency-prediction  programme  (referred  to  in  Section  3.2). 

It  was  found  that  a reduction  of  the  second  stage  flow  areas  made  the  second  rotor  incidence 
angles  more  positive  but  the  first  rotor  incidence  angles  more  negative,  with  a small  net  loss 
in  mean  total-total  efficiency,  and  a marked  decrease  in  mean  total-static  efficiency  because  of 
the  increased  exit  kinetic  energy  loss.  Conversely,  it  was  found  that  an  increase  of  the  second 
stage  flow  areas  decreased  the  mean  total-total  efficiency  but  increased  the  mean  total-static 
efficiency  as  a result  of  the  reduced  exit  kinetic  energy  loss  associated  with  reduced  stage  loading. 
Both  results  are  clearly  indicated  in  Fig.  12. 

The  above  attempt  to  obtain  improved  mean  efficiencies  was  inconclusive  since  changes  in 
first  rotor  incidence  angles  and  second  stage  exit  kinetic  energy  masked  any  beneficial  effect 
of  incidence  angle  change  on  the  second  rotor:  only  by  specifically  varying  the  second  rotor 
inlet  blade  angles  could  the  true  effect  be  estimated. 

5.3  Comparison  of  Estimated  Mean  Efficiencies  of  Axial  Turbines  with  Estimated  Mean 
Efficiency  of  Radial  Turbine 

The  present  study  with  axial-flow  turbines  is  an  extension  of  the  investigation  performed 
on  the  radial-inflow  turbine  (Erm  (1973))  to  determine  the  best  turbine  configuration  for  operation 
with  the  pulsating  flow  from  an  experimental  combustor.  Consequently,  it  is  appropriate  to 
compare  the  results  obtained  for  the  two  types  of  turbine. 

The  estimated  mean  efficiencies  of  the  two-stage  axial-flow  turbines  considered  in  this 
investigation  are  listed  in  Section  5.2.3.  Build  3,  which  has  a maximum  mean  total-static  efficiency 
of  about  75  -6%,  gives  the  best  results.  The  corresponding  maximum  mean  total-static  efficiency 
of  the  radial-inflow  turbine  is  about  76  -3%.  Thus  the  two  types  of  turbine  have  much  the  same 
efficiency.  It  is  to  be  noted,  however,  that  a number  of  other  factors  have  to  be  considered 
before  any  firm  conclusions  can  be  drawn  on  the  relative  merits  of  the  different  turbines. 

Although  the  general  form  of  the  pulsating  flow  used  in  both  investigations  is  the  same 
(blowdown  with  uniform  combustor  conditions  and  scavenging  in  which  the  incoming  scavenging 
air  intermixed  completely  with  the  undischarged  combustor  contents),  and  although  an  attempt 
was  made  to  maintain  a comparable  severity  of  pulsation  (refer  to  Section  4),  the  two  types 
of  pulsating  flow  are  in  actual  fact  different  due  to  the  fact  that  the  temperatures  and  pressures 
at  the  important  cycle  locations  are  not  the  same  in  both  cases. 

A factor  that  must  be  considered  when  comparing  the  performance  of  the  axial-flow  turbines 
with  that  of  the  radial-flow  type  is  the  relative  sizes  of  the  two  different  types  of  machine.  The 
capacities  of  the  axial-flow  types  are  smaller  than  that  of  the  radial-flow  type,  so  the  axial-flow 
turbines  may  have  in  fact  suffered  some  penalty  as  a consequence. 

Another  factor  that  must  be  considered  when  comparing  the  performance  of  the  two  types 
of  turbine  is  that  no  diffuser  was  incorporated  in  the  axial-flow  turbine  study  whereas  one  was 
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used  with  the  radial  machine.  Consequently,  absolute  values  of  mean  efficiency  for  the  two 
types  are  not  directly  comparable. 

The  choice  between  a two-stage  axial-flow  turbine  and  a single-stage  radial-inflow  turbine 
will  clearly  not  be  decided  on  turbine  efficiency  alone.  Since  the  performance  is  much  the  same, 
it  is  likely  that  reliability,  maintenance  needs  and  cost  may  be  more  decisive  factors;  these 
factors  are  all  generally  improved  by  simplicity  of  construction. 

6.  CONCLUSIONS 

An  endeavour  has  been  made  to  estimate  the  mean  efficiencies  of  a number  of  different 
two-stage  axial-flow  turbines  under  pulsating-flow  operating  conditions.  Initially  the  turbine 
loss  characteristics  were  determined  by  using  the  revised  Ainley  and  Mathieson  correlation  but 
it  was  soon  found  that  this  correlation  was  inadequate  at  high  negative  angles  of  incidence  and 
consequently  some  modification  was  needed  to  the  correlation  in  this  area  of  operation.  It 
became  quite  clear  that  a better  knowledge  of  turbine  loss  characteristics  at  high  negative  angles 
of  incidence  will  be  needed  before  more  precise  optimisation  of  turbine  geometry  for  operation 
under  pulsating  flow  can  be  attempted. 

The  maximum  mean  total-static  efficiencies  that  could  reasonably  be  expected  in  practice 
from  the  four  basic  turbine  configurations  under  pulsating-flow  operation  are  listed  below. 


Build  1 

74-7% 

Build  2 

70-1% 

Build  3 

75-6% 

Build  4 

70-7% 

A turbine  with  outlet  blade  angles  in  the  range  50°  to  60°  in  the  first  stage  and  outlet  blade 
angles  in  the  range  60°  to  70"  in  the  second  stage  could  possibly  give  better  mean  total-static 
efficiencies  under  pulsating-flow  conditions  than  the  four  turbines  investigated. 

An  analysis  of  the  flow  through  a turbine  during  a pulse  has  shown  that  the  mean  inlet 
flow  angles  encountered  by  the  different  blade  rows  varies  and  it  may  be  possible  to  obtain 
improved  values  of  mean  efficiency  if  the  matching  between  the  inlet  flow  angles  and  the  inlet 
blade  angles  throughout  the  turbines  is  improved.  An  attempt  to  show  this  was  inconclusive. 

Since  the  indicated  performances  of  a single-stage  radial-inflow  turbine  and  a two-stage 
axial-flow  turbine  for  operation  with  the  assumed  pulsating-flow  conditions  are  much  the  same, 
then  the  choice  between  these  two  types  of  turbine  will  be  governed  by  factors  such  as  simplicity 
of  construction  and  cost  of  the  engine  rather  than  by  performance. 
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APPENDIX  I 


r 


PRINCIPAL  NOTATION 

A flow  area  measured  normal  to  camber  line,  m* 

A a annulus  area,  m* 

Ay  flow  area  measured  normal  to  flow  direction,  m* 

a acoustic  velocity,  m/s 

c blade  chord,  m 

Cp  specific  heat  at  constant  pressure.  J/kg.K 

h annulus  height,  m 

i incidence  angle  of  flow  on  to  a blade  row  given  by  difference  between  gas  flow 

angle  relative  to  blade  inlet  and  blade  inlet  angle,  rad 

U stalling  incidence  angle,  rad 

k radial  tip  clearance,  m 

M Mach  number. 

mu  mass  of  gas  in  combustor,  kg 

m mass  flow  rate,  kg/s 

nioi  t mass  flow  rate  of  gas  out  of  combustor,  kg/s 
A/min  increment  of  mass  flowing  into  combustor,  kg 

Amot  t increment  of  mass  flowing  out  of  combustor,  kg 

P total  or  stagnation  pressure.  Pa 

p static  pressure.  Pa 

Pscav  presawe  of  scavenging  air  at  entry  to  combustor.  Pa 

PR tit7  turbine  or  stage  pressure  ratio  based  upon  total  pressure  at- station -1-and-total 

pressure  at  station  7. 

P Rt7sis  stage  pressure  ratio  based  upon  total  pressure  at  station  7 and  static  pressure 
at  station  13. 

P Rt7tis  stage  pressure  ratio  based  upon  total  pressure  at  station  7 and  total  pressure 
at  station  13. 

P Rtisi3  turbine  pressure  ratio  based  upon  total  pressure  at  station  1 and  static  pressure 
at  station  13. 

R gas  constant.  J/kg.K 

s blade  pitch,  m 

T total  or  stagnation  temperature.  K 

t static  temperature.  K 

Tscav  temperature  of  scavenging  air  at  entry  to  combustor.  K 
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fuMAssi  mean  combustor  cycle  temperature  computed  on  a mass-averaged  basis.  K 
fuTiME)  mean  combustor  cycle  temperature  computed  on  a time-averaged  basis.  K. 
te  blade  trailing  edge  thickness  measured  normal  to  camber  line  at  trailing  edge,  m 

rmkX  maximum  blade  thickness,  m 


A / time  increment,  s 

U rotor  blade  speed  at  mean  diameter,  m/s 

V gas  velocity,  m/s 

^stis7  spouting  velocity  based  upon  total  pressure  at  station  1 and  static  pressure  at 
station  7.  m/s 

v0  volume  of  combustor.  ms 

W specific  shaft  work.  J/kg 

A If  act  actual  increment  of  work  resulting  from  flow  of  an  increment  of  mass  through 

turbine.  J 

Alf'is  isentropic  or  theoretical  increment  of  work  resulting  from  flow  of  an  increment 
of  mass  through  turbine.  J 

Yr  profile  loss  coefficient. 

Ya  secondary  loss  coefficient. 

Yk  tip  clearance  loss  coefficient. 

Kt  total  loss  coefficient. 


Greek  Letters 

a gas  flow  angle  measured  relative  to  axial  direction,  rad 

fi  blade  angle  measured  relative  to  axial  direction,  rad 

Y specific  heat  ratio. 

riTisv  turbine  or  stage  efficiency  based  upon  total  pressure  at  station  1 and  static 
pressure  at  station  7. 

7JT1T7  turbine  or  stage  efficiency  based  upon  total  pressure  at  station  1 and  total 
pressure  at  station  7. 

71T7T13  stage  efficiency  based  upon  total  pressure  at  station  7 and  total  pressure  at 
station  13. 

t)tisi3  turbine  efficiency  based  upon  total  pressure  at  station  1 and  static  pressure 
at  station  13. 

t)titi3  turbine  efficiency  based  upon  total  pressure  at  station  I and  total  pressure  at 
station  13. 

tjtisis  mean  turbine  efficiency  for  cyclic  pulsating-flow  operating  conditions  based 
upon  total  pressure  at  station  1 and  static  pressure  at  station  1 3. 

7JTIT13  mean  turbine  efficiency  for  cyclic  pulsating-flow  operating  conditions  based 
upon  total  pressure  at  station  I and  total  pressure  at  station  1 3. 

vns7  velocity  ratio  based  upon  total  pressure  at  station  1 and  static  pressure  at  station  7. 

vTisis  velocity  ratio  based  upon  total  pressure  at  station  1 and  static  pressure  at 
station  13. 
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VT1T1S 


velocity  ratio  based  upon  total  pressure  at  station  1 and  total  pressure  at 
station  13. 

stage  loading  coefficient  ('F  = AW/U1). 


Subscripts 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 
13 

V,  W 


station  within  the  combustor. 

station  at  first  stator  inlet. 

station  at  first  stator  throat. 

station  downstream  of  first  stator. 

station  at  first  rotor  inlet. 

station  at  first  rotor  throat. 

station  downstream  of  first  rotor. 

station  at  second  stator  inlet. 

station  at  second  stator  throat. 

station  downstream  of  second  stator. 

station  at  second  rotor  inlet. 

station  at  second  rotor  throat. 

station  downstream  of  second  rotor. 

station  downstream  of  turbine. 

locations  on  pulsed  combustor  cycle  (see  Fig.  8). 
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APPENDIX  II 


A2.  CALCULATION  OF  TURBINE  LOSS  CHARACTERISTICS 

The  technique  used  to  calculate  the  turbine  loss  characteristics  by  means  of  the  revised 
Ainley  correlation  is  now  described.  The  analysis  is  presented  for  blade  row  456  of  the  turbine 
under  investigation  (refer  to  Fig.  I)  but  is  applicable  to  any  arbitrary  blade  row  provided  that 
the  correct  geometrical  dimensions  are  used  and  the  subscripts  are  modified  appropriately. 
It  is  also  to  be  noted  that  the  following  analysis  is  presented  for  subsonic  flow  conditions  except 
where  otherwise  indicated. 

A2. 1 Profile  Losses 

The  profile  loss  coefficient  at  zero  incidence  can  be  determined  by  means  of  the  relationship 


-[ 


Yp4«0)  = I lp<O)(0A-O)  + 


m 


Fp(0)  <o- — <*.,  — Fp(0) 


(A2.1) 


In  this  equation,  Fp(0(  (e4-  oi  is  the  profile  loss  coefficient  at  zeio  incidence  for  a blade 
having  (34  = 0°  and  same  as  and  s/c  as  the  actual  blade.  Kp(o>  ((i4-  -«6,  is  the  profile  loss  coeffi- 
cient at  zero  incidence  for  a blade  having  P«  = — a5  and  same  a6  and  s/c  as  the  actual  blade. 
The  values  of  Fp(0)  (o4-  oi  and  Fp(oi  <o4~  -»6i  can  be  determined  from  graphs  given  by  Ainley 
and  Mathieson  (1951). 

In  order  to  determine  the  profile  loss  coefficient  for  a particular  angle  of  incidence,  i,  it  is 
first  necessary  to  determine  the  stalling  incidence  angle,  /»,  from  graphs  given  by  Ainley  and 
Mathieson  (1951)  and  then,  after  computing  i/ia , determine  yp4«/Fp4#«»  from  another  graph 
given  in  the  same  reference.  Once  Ynel  Ff48<o>  has  been  determined  then  it  is  a simple  matter 
to  calculate  Fp48. 

For  supersonic  flow  conditions  only,  it  is  necessary  to  use  an  additional  term, 
viz.  [1  -t- 60 (Mt—  1)*],  when  determining  profile  losses.  This  term  has  not  been  used  in  the 
current  investigation  since  all  flow  conditions  considered  have  been  subsonic. 

A2.2  Secondary  Losses 

The  secondary  loss  coefficient  for  any  specified  angle  of  incidence  can  be  determined  by 
means  of  the  relationship 

Fs4«  =0  0334 
Ci 


where 

and 


s/c 


XCOS  a6l  r Ci.'l2T  COS2  1 

cos  P4J  [_.r/cj  LcossaM46J 

= 2 (tan  a4  — tan  a5)  cos  aM46 


tan  aM46  = 0-5  (tan  a4  + tan  a5) 


(A2.2) 

(A2.3) 

(A2.4) 


A2.3  Tip  Clearance  Losses 

The  tip  clearance  loss  coefficient  for  any  specified  angle  of  incidence  can  be  determined 
by  means  of  the  relationship 


Y K46  = B 


vi  m ° tcl]2  r cos2  «s  I 

aIW  Lv/cJ  Lcos,*m«J 


(A2.5) 
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The  value  of  B is  0-47  for  unshrouded  blades  while — and  aM46  can  be  determined  from 

i/c 

equations  (A2.3)  and  (A2.4)  respectively. 


A2.4  Total  Losses 

The  total  loss  coefficient  for  any  specified  angle  of  incidence  can  be  determined  by  simply 
adding  together  the  component  loss  coefficients  for  this  angle  of  incidence  as  follows: 

Tt46  = y P46  + Ts4«  + y K46  (A2.6) 


The  variation  of  the  component  loss  coefficients,  Tp46,  Ts4«  and  Ykm,  as  well  as  the  total 
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A3.  THERMODYNAMIC  ANALYSIS  OF  FLOW  THROUGH  TURBINE  AND 

PREDICTION  OF  TURBINE  EFFICIENCY 

The  details  of  the  equations  used  and  the  technique  adopted  in  order  to  establish  the  flow 
conditions  through  the  turbine  to  enable  the  efficiency  to  be  predicted  are  now  described. 

A3.1  Known  Data  Used  in  Analysis 

The  known  data  used  in  the  analysis  are  listed  in  a number  of  different  groups  as  follows. 

A3. 1.1  Specified  Turbine  Operating  Conditions 

The  specified  turbine  operating  conditions 

Tx 
Px 
Pi 
U 

A3. 1.2  Known  Thermodynamic  Data 

The  known  thermodynamic  data  are  listed  below.  It  is  assumed  that  the  specific  heats  and 
other  gas  constants  do  not  vary  with  temperature. 

Cp  =1-00  kJ/kg.K 
R = 0-287  kJ/kg.K 

Y = 1-4 

A3. 1 .3  Known  Turbine  Geometrical  Data 

The  known  turbine  geometrical  data  are  listed  below. 

A2  — 1529  mm2 

/fA3  = 3057  mm2 

Ab  = 1822  mm2 

Aa,  — 3644  mm2 

P2  = —60°  (Nominal) 

p5  = —60°  (Nominal) 

It  should  be  noted  that  Ainley  and  Mathieson  (1951)  give  a method  of  calculating  the  gas 
outlet  angle  as  a function  of  blade  geometrical  characteristics  and  flow  Mach  number.  Since 
it  is  shown  that  for  straight-backed  blades  and  M > 0-5  the  actual  gas  outlet  angle  does  not 
vary  more  than  a few  degrees,  it  was  decided  to  ignore  this  variation  for  the  purpose  of  the 
present  study.  Thus  a blade  design  is  assumed  which  will  give  the  nominated  gas  outlet  angle 
at  the  higher  Mach  numbers,  and  no  variation  of  gas  outlet  angle  with  conditions  is  allowed 
for.  Revision  of  the  method  to  include  gas  outlet  angle  deviation  with  flow  conditions  could  be 
incorporated  for  greater  accuracy,  provided  the  accuracy  of  other  data  (loss  coefficients,  etc.) 
is  good  enough  to  warrant  this  improvement. 

A3. 2 Details  of  Equations  Used  to  Establish  Flow  Conditions  Through  Turbine 

As  has  been  stated  in  Section  2.2,  in  order  to  establish  the  flow  conditions  through  the 
turbine  for  the  chosen  operating  conditions  it  is  necessary  to  adopt  an  iterative  approach.  This 
process  involves  systematically  choosing  a series  of  trial  values  of  ps  and  evaluating  tentative 


are  listed  below. 
K 
Pa 
Pa 
m/s 


values  of  the  flow  variables  through  the  turbine  until  the  mass  flow  rates  through  the  stator 
and  rotor  agree  within  defined  limits. 

It  is  necessary  to  use  alternative  sets  of  equations  to  evaluate  the  flow  variables  through 
the  turbine  depending  upon  whether  or  not  the  flow  is  subsonic  or  supersonic.  Both  sets  of 
equations  are  considered  in  the  following. 


A3. 2.1  Analysis  of  Flow  Through  Stator 

The  analysis  of  the  flow  through  the  stator  is  subdivided  into  a number  of  different  sections 
each  of  which  is  an  analysis  of  the  flow  at  a particular  stator  station  for  either  subsonic  or 
supersonic  flow  conditions  downstream  of  the  stator. 


A3. 2. 1.1  Subsonic,  sonic  or  supersonic  flow  downstream  of  stator  (any  \f3) 

The  following  equations,  viz.  (A3.1)  to  (A3.7),  used  to  determine  flow  variables  are  applicable 
for  any  type  of  flow  downstream  of  the  stator. 

Stat  'on3:  Downstream  of  stator 

The  total  temperature  is  given  by  the  following  relationship: 

T3  = Tx  (A3.1) 


The  total  loss  coefficient  for  the  stator  is 

Ttis=  (A3.2) 

P*  ~P* 

This  relationship  can  be  transposed  to  yield  the  following  expression  for  the  total  pressure: 

r.  - V/™'’1  (A3.3) 

1 + JT13 

Now  that  T3,  P3  and  p3  are  all  known,  it  is  possible  to  determine  the  static  temperature 
by  means  of  the  isentropic  relationship 


—ir 

It  is  now  also  possible  to  determine  the  absolute  gas  velocity  as  follows: 


(A3.4) 


while  the  acoustic  velocity  is 

y3  = V[2C„(T3  - /,)] 

(A3. 5) 

so  yielding  the  Mach  number 

a3  = VM'j] 

y. 

(A3.6) 

A f,  =- 

(A3. 7) 

If  M3<  10  then  the  flow  at  the  throat  is  subsonic  while  if  M3  > 1 0 then  the  flow  at 
the  throat  is  sonic.  The  equations  applicable  to  both  of  these  flow  conditions  are  considered 
in  turn. 

A3. 2.1. 2 Subsonic  flow  downstream  of  stator  (A/,  < 1 -0) 


The  following  equations,  viz.  (A3. 8)  to  (A3. 18),  used  to  determine  flow  variables  are 
applicable  when  the  flow  downstream  of  the  stator  is  subsonic.  For  this  flow,  conditions  at 
stations  2 and  3 are  essentially  the  same. 

Station  2:  Stator  throat 


Flow  variables  at  the  stator  throat  are  given  by  the  following  relationships: 


= r, 

(A3. 8) 

= '. 

(A3.9) 

= p* 

(A3. 10) 

= p* 

(A3. 11) 

= y. 

(A3. 12) 

= «« 

(A3. 13) 
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m,  = m3 
<**  = Pi 

The  mass  flow  rate  can  be  determined  as  follows: 

^ PtAtVt 


(A3. 14) 
(A3. 15) 


(A3. 16) 


Station  3 : Downstream  of  stator 

Flow  variables  downstream  of  the  stator  are  given  by  the  following  relationships: 

(A3. 17) 

“>  = *«  (A3. 18) 

A3.2.1.3  Sonic  or  supersonic  flow  downstream  of  stator  (A/,  ^ 1 -0) 

The  following  equations,  viz.  (A3. 19)  to  (A3.29),  used  to  determine  flow  variables  are 
applicable  when  the  flow  downstream  of  the  stator  is  sonic  or  supersonic.  For  supersonic  flow, 
conditions  at  stations  2 and  3 are  different  as  a consequence  of  a Prandtl-Meyer  expansion 
between  these  stations. 

Station  2 : Stator  throat 

Flow  variables  at  the  stator  throat  are  given  by  the  following  relationships: 

T%  = Ti  (A3. 19) 

pi  = p*  (A3. 20) 

= P*  (A3. 21) 

A relationship  between  the  total  temperature  and  static  temperature  is 

T r.  . (r  — 1) 1 


(A3. 22) 


Since  the  velocity  at  the  stator  throat  is  sonic  = 1 -0),  then  the  equation  for  the  static 
temperature  is 

(A3a) 


A corresponding  equation  for  the  static  pressure  is 


b^Y 


The  absolute  gas  velocity  can  be  determined  from  the  relationship 

V*  = \/[2Cp(rl  - rt)] 

The  mass  flow  rate  can  now  be  determined  from 

■ _ Pt  a,  y2 
2 rtf2 

Station  3:  Downstream  of  stator 

The  mass  flow  rate  is  given  by  the  following  relationship: 


Now,  the  mass  flow  rate  can  be  related  to  flow  conditions  as  follows: 

. _PiA„  F,cosa, 


This  expression  can  be  transposed  to  yield 

*»  — ± arcos 


y.p*Aktv%\ 


(A3. 24) 


(A3. 25) 


(A3. 26) 


(A3. 27) 


(A3. 28) 


(A3. 29) 


where  t3  and  Va  are  given  by  equations  (A3.4)  and  (A3.5)  respectively. 

In  accordance  with  the  sign  convention  adopted  for  gas  flow  angles,  the  sign  of  a,  can  be 
determined  as  follows: 

if «,  > 0 then  as  > 0 
if  a,  < 0 then  a,  ^ 0 


A3.2.2  Analysis  of  Flow  at  Rotor  Inlet 

The  following  equations,  viz.  (A3.30)  to  (A3.40),  used  to  determine  flow  variables  at  the 
rotor  inlet  are  applicable  regardless  of  Mach  number  downstream  of  the  stator.  The  subscripts  4, 
S and  6 on  flow  variables  indicate  that  these  variables  are  specified  relative  to  the  rotor. 

Flow  variables  at  the  rotor  inlet  are  given  by  the  following  relationships: 


'4  = '» 
Pi  =Ps 
m«  = m, 
a«  = a3 


With  reference  to  Fig.  1 it  is  apparent  that 

V3  cos  a,  = Vt  cos  a4 

and,  having  regard  for  the  sign  convention, 

V3  sin  a,—  U — — V3  sin  a4 

so  therefore 


and 


*4  = 


— artan 


V3  sin  a,  — t/~ 
y3  cos  otj 


y3  cos  a, 
cos  a4 


The  total  temperature  can  be  determined  from  the  relationship 


The  total  pressure  is  given  by  the  isentropic  relationship 


The  Mach  number  can  be  determined  as  follows: 


Mi 


(A3. 30) 
(A3. 31) 
(A3.32) 
(A3. 33) 


(A3. 34) 


(A3. 35) 


(A3. 36) 


(A3. 37) 


(A3. 38) 


(A3. 39) 


(A3.40) 


A3. 2. 3 Analysis  of  Flow  Through  Rotor 

The  analysis  of  the  flow  through  the  rotor  is  subdivided  into  a number  of  different  sections 
each  of  which  is  an  analysis  of  the  flow  at  a particular  rotor  station  for  either  subsonic  or  super- 
sonic flow  conditions  downstream  of  the  rotor. 

Since  the  analysis  of  the  flow  through  the  rotor  is  very  similar  to  the  analysis  considered 
earlier  for  the  stator,  in  order  to  simplify  the  presentation,  the  equations  used  to  determine 
flow  variables  will  often  just  be  listed  without  the  accompanying  text  and  intermediate  steps. 


A3. 2. 3.1  Subsonic,  sonic  or  supersonic  flow  downstream  of  rotor  ( any  A/,) 

The  following  equations,  viz.  (A3.41)  to  (A3.47),  used  to  determine  flow  variables  are 
applicable  for  any  type  of  flow  downstream  of  the  rotor. 
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Station  6:  Downstream  of  rotor 

Flow  variables  downstream  of  the  rotor  are  given  by  the  following  relationships: 

T,  = Tt 
P»  =Pi 

where  />,  »«  a specified  turbine  operating  condition. 

Ot  ow  variables  downstream  of  the  rotor  can  be  determined  as  follows: 

_ _ Pt  + YnuPt 
* 1 + yT4« 

where  It4«  is  a function  of  the  incidence  angle,  »4 


(A3.41) 
(A3 .42) 


(A3. 43) 


(A3.44) 

y.  = V[2Cp(r,  - r4)] 

(A3.45) 

ai  = v/[r*'«] 

(A3. 46) 

(A3.47) 

If  Mt  < 10  then  the  flow  at  the  throat  is  subsonic  while  if  Mt  > 10  then  it  is  sonic. 
The  equations  applicable  to  both  of  these  flow  conditions  shall  be  considered  in  turn. 

A3. 2.3.2  Subsonic  flow  downstream  of  rotor  (A/*  < 1 -0) 

The  following  equations,  viz.  (A3.48)  to  (A3. 58),  used  to  determine  flow  variables  are 
applicable  when  the  flow  downstream  of  the  rotor  is  subsonic.  For  this  flow,  conditions  at 
stations  5 and  6 are  essentially  the  same. 

Station  5 : Rotor  throat 

Flow  variables  at  the  rotor  throat  are  given  by  the  following  relationships: 

r,  = r4  (A3.48) 

t>  = /.  (A3.49) 

Pi  = P»  (A3. 50) 

/>»=/>.  (A3. 51) 

Vs  = K,  (A3.52) 

at  = a,  (A3.53) 

Mi  = M,  (A3. 54) 

a,  = p,  (A3. 55) 

The  mass  flow  rate  can  be  determined  by  means  of  the  relationship 

(A3-56> 

Station  6:  Downstream  of  rotor 

Flow  variables  downstream  of  the  rotor  are  given  by  the  following  relationships: 

mt  — m,  (A3.57) 

*4  = at  (A3.58) 

A3.2.3.3  Sonic  and  supersonic  flow  downstream  of  rotor  (Mt  > 1 0) 

The  following  equations,  viz.  (A3.59)  to  (A3.67),  used  to  determine  flow  variables  are 
applicable  when  the  flow  downstream  of  the  rotor  is  sonic  or  supersonic. 


Station  5:  Rotor  throat 

Flow  variables  at  the  rotor  throat  are  given  by  the  following  relationships : 

Tt  = Tt  (A3. 59) 

Pi  = Pi  (A3.60) 

= P,  (A3.6I) 

Other  flow  variables  at  the  rotor  throat  can  be  determined  by  means  of  the  relationships 
shown  below : 


/ r» 

W] 


Pi  = 


FW' 


Vi  = Vl2Cp(r5 

Pi  A i F, 


m.  = 


ft's 


Station  6:  Downstream  of  rotor 

The  mass  flow  rate  is  given  by  the  following  relationship: 

mt  = ni. 

The  gas  flow  angle  can  be  determined  by  means  of  the  relationship  shown  below: 


«»  = ± arcos 


ft '«! 
_P*A  as*VJ 


(A3.62) 

(A3.63) 

(A3.64) 
(A3. 65) 

(A3. 66) 

(A3. 67) 


where  f,  and  Vt  are  given  by  equations  (A3.44)  and  (A3.45)  respectively. 

In  accordance  with  the  sign  convention  adopted  for  gas  flow  angles,  the  sign  of  a6  can  be 
determined  as  follows: 

if  as  > 0 then  a,  > 0 
if  as  $ 0 then  a,  < 0 

A3.2.4  Analysis  of  Flow  at  a Fixed  Station  Downstream  of  Rotor 

Since  the  analysis  of  the  flow  at  a fixed  station  downstream  of  the  rotor  is  very  similar 
to  the  analysis  considered  earlier  for  the  rotor  inlet,  then  in  order  to  simplify  the  presentation, 
the  equations  used  to  determine  flow  variables  will  often  just  be  listed  without  the  accompanying 
text  and  intermediate  steps. 

The  following  equations,  viz.  (A3.68)  to  (A3.75),  used  to  determine  flow  variables  at  a fixed 
station  downstream  of  the  rotor  are  applicable  regardless  of  Mach  number  downstream  of  the 
rotor. 


V , 


(A3.68) 

= m. 

(A3. 69) 

= «• 

(A3. 70) 

f F.sin  a,-  U~\ 

artan[  K,cos«,  J 

(A3. 71) 

cos  *, 

?cos«7 

J 

(A3. 72) 

T IP  t 

(A3.73) 
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r 


T 


(A3.74) 
(A3. 75) 


t 


i 


A3.3  Balancing  of  Mass  Flow  Rates 

As  has  been  indicated  in  Section  A3.2  in  order  to  establish  the  flow  conditions  through 
the  turbine  for  the  chosen  operating  conditions  it  is  necessary  to  adopt  an  iterative  approach. 
This  process  involves  systematically  choosing  a series  of  trial  values  of  pt  and  evaluating  the 
flow  variables  through  the  turbine  until  the  mass  flow  rates  through  the  stator  and  rotor,  as 
determined  by  equations  (A3.16)  or  (A3.26)  (stator)  and  (A3.56)  or  (A3.65)  (rotor),  agree  within 
defined  limits. 


A3. 4 Prediction  of  Turbine  Efficiency 

As  soon  as  the  flow  conditions  through  the  turbine  have  been  firmly  established  for  the  chosen 
operating  conditions  it  is  possible  to  determine  the  total-static  efficiency  and  the  velocity  ratio 
as  well  as  the  specific  work. 

The  total-static  efficiency  can  be  determined  as  follows: 


7JT1S7  = 


The  velocity  ratio  is  given  by 


(A3. 76) 


VT1S7  = 


U 

FST1S7 


(A3. 77) 


where  Fstis7  is  the  spouting  velocity  which  is  defined  as  that  velocity  which  has  an  associated 
kinetic  energy  equal  to  the  isentropic  enthalpy  drop  across  the  turbine.  If  Kstis?  is  expressed 
in  terms  of  thermodynamic  variables,  the  velocity  ratio  becomes 


VT1S7 


U 


Finally,  the  specific  work  can  be  determined  as  follows : 

^17  = Q (Ty  - T,) 


(A3. 78) 


(A3. 79) 
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A4.  ESTIMATION  OF  TURBINE  MEAN  EFFICIENCY  UNDER  PULSATING-FLOW 

CONDITIONS 

The  technique  used  to  estimate  the  mean  efficiencies  of  the  two-stage  turbines  when  they 
are  operating  under  the  idealised  pulsating-flow  conditions  outlined  in  Section  4 and  Fig.  8 
is  now  described.  In  the  following  analysis  static  pressures  and  temperatures  are  used  when 
dealing  with  combustor  conditions.  It  is  assumed  that  total  pressures  and  temperatures  at  the 
turbine  inlet  are  equal  to  their  respective  static  values  within  the  combustor.  The  analysis  is 
simplified  by  assuming  that  all  flow  is  one-dimensional  and  that  the  specific  heats  and  other 
gas  constants  do  not  vary  with  temperature. 

A4. 1 Blowdown  Phase 

During  the  blowdown  phase  there  is  no  cold  scavenging  air  flowing  into  the  combustor 
and  the  pressure  continually  decreases  as  the  hot  gas  in  discharged  through  the  turbine. 

The  analysis  of  the  blowdown  phase  is  described  for  the  first  time  interval  (01),  but  is 
applicable  to  all  subsequent  intervals. 

In  order  to  commence  the  calculations  it  is  necessary  to  select  a trial  value  of  fov,  the 
combustor  temperature  before  combustion  (see  Fig.  8). 

The  initial  mass  of  gas  in  the  combustor  is  given  by 


'Wo(O)  = (A4.1) 

i\  fov 

Instead  of  computing  the  true  mass  flow  rates  through  the  different  turbines,  a simpler 
approach  was  adopted  whereby  pseudo  mass  flow  rates  were  determined  after  first  assuming  that 
the  turbines  were  replaced  by  a nozzle,  of  area  At,  acted  upon  by  the  pressure  ratio  Pi(0)/pi».  With 
this  simplified  approach  the  relationship  between  the  mass  discharged  from  the  combustor,  the 
combustor  temperature  and  pressure  (i.e.  turbine  inlet  conditions)  and  the  turbine  efficiency,  at 
any  instant  of  the  pulsating-flow  cycle,  is  the  same  as  if  the  true  mass  flow  rates  had  been  used 
and,  consequently,  the  computed  mean  efficiencies  are  the  same  as  those  which  would  be  obtained 
by  using  the  true  mass  flow  rates.  For  the  simplified  approach,  the  critical  pressure  ratio,  i.e.  the 
pressure  ratio  corresponding  to  Mt(0)  = 10,  can  be  determined  as  follows: 


{'ll  =ri+^i 

(PiaJcRIT  L 2 


If 


If 


{- 

(.Pm 


P.(  0) 

Pis 

Mt{  0)  = 1 0 
P,(0) 

Pis 


1 +• 

i}  CRIT 


Y 

Y-t 


o 

- .1 


then 


then 


CRIT 


(Y-l) 


(A4.2) 


(A4.3) 


(A4.4) 


Now  that  the  value  of  Af,(0)  is  known,  the  pseudo  mass  flow  rate  through  the  turbine  (out 
of  the  combustor)  can  be  determined  as  follows : 
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moiMO)  = AiM^0)Pl(0) 


(A4.5) 


/ Y 

r i i 

/ rtm 

if 

5- 

1 <N 

+ 

1 

Y+l 

«Y-1) 


This  rate  is  assumed  to  apply  for  the  whole  interval,  and  the  increment  of  outflow  for 
the  period  is 


A/Woi>t(01)  = mouT(0)-Af(01) 


(A4.6) 


The  ideal  or  isentropic  increment  of  work  obtainable  from  the  turbine  in  the  interval  is 
therefore 


A^is(Ol)  = AmourCOO.Cp.r^O) 


(A4.7) 


In  order  to  determine  the  actual  increment  of  work,  the  initial  turbine  efficiency,  ijTisisfO), 
must  be  determined.  Before  this  can  be  done,  it  is  necessary  to  calculate  the  initial  velocity  ratio, 
vtisis(0),  by  means  of  the  following  relationship: 


vtisis(0)  = 


U 


2CP7\(0) 


L >i(o)j 


!]]] 


(A4.8) 


where  5 is  the  number  of  stages. 

The  value  of  y)Tisis(0)  is  determined  from  the  appropriate  y]tisi3  versus  vTisis  relationship 
shown  in  Fig.  7.  It  is  now  possible  to  determine  the  actual  increment  of  work  as  follows: 


AH/act(01)=  A W/is(01).y)T1si3(0) 


(A4.9) 


Before  the  flow  through  the  turbine  during  time  interval  (12)  can  be  analysed,  the  combustor 
pressure  and  temperature  at  the  end  of  interval  (01)  must  be  determined.  The  mass  of  gas  in 
the  combustor  at  the  end  of  interval  (01)  is 


m„(l)  = m„(0)  — Amou-KOl) 


(A4.10) 


Assuming  that  the  gas  remaining  in  the  combustor  expands  isentropically  during  interval 
(01),  the  new  conditions  of  pressure  and  temperature  are 

(A4") 

These  are  the  initial  conditions  for  time  interval  (12). 

The  flow  through  the  turbine  during  time  interval  (12)  can  now  be  analysed.  The  procedure 
used  is  the  same  as  that  outlined  above  for  time  interval  (01)  except  that  T,(0)  is  replaced  by 
T^l),  and  likewise  for  other  variables.  This  process  can  be  repeated  for  subsequent  time  intervals 
until  the  combustor  pressure  is  found  to  be  equal  to  or  less  than  the  scavenging  pressure,  Fscav- 
When  this  occurs  the  scavenging  phase  is  commenced. 


A4.2  Scavenging  Phase 

During  the  scavenging  phase  the  cold  incoming  air  is  assumed  to  intermix  completely  with 
the  hot  undischarged  combustor  contents  and  the  pressure  in  the  combustor  remains  constant. 

For  the  first  time  interval  (//)  of  the  scavenging  phase,  the  turbine  entry  conditions,  7\(/) 
and  /*,(/),  are  equal  to  those  at  the  end  of  blowdown  and  therefore  Amoi:i(/J),  AIFis(/7)and 
SWinOJ)  can  be  determined  by  using  the  same  method  as  that  outlined  above  for  the  blow- 
down phase.  During  the  first  time  interval  (A/)  the  increment  of  mass  of  cold  scavenging  air 
flowing  into  the  combustor  is  given  by  the  following  relationship  (volume  inflow  equals  volume 
outflow): 
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a-  l^I)  

hmiN(IJ)  — — ^ (A4. 1 3) 

Before  the  flow  through  the  turbine  during  time  interval  ( JK ) can  be  analysed,  it  is  necessary 
to  determine  the  combustor  temperature  at  the  end  of  interval  (//).  An  expression  for  this 
temperature  is  now  developed. 

At  the  end  of  time  interval  (IJ)  the  combustor  contains  [m0(I)  — Aw0l  t(A/)]  kg  of  hot  gas 
at  temperature  /„(/)  and  \mm(U)  kg  of  cold  scavenging  air  at  temperature  Tscav-  That  is 

m«(J)  = m^I)  - Amo u !</./)  + hmiV(IJ)  (A4.14) 

If  the  hot  gas  and  the  cold  air  mix  completely,  then  the  temperature  of  the  mixture  is  given 
by  the  relationship: 


to U)  = 


/of/)  [m0(7)  — Amorr(/-/)]  + TscAV.&ntwifJ) 


(A4.15) 


This  is  the  initial  temperature  for  time  interval  (JK). 

The  flow  through  the  turbine  during  time  interval  (JK)  can  now  be  analysed.  The  procedure 
used  is  the  same  at  that  outlined  above  for  time  interval  (IJ)  except  that  T,(/)  is  replaced  by 
Ti(J)  and  likewise  for  other  variables.  This  process  can  be  repeated  for  subsequent  time  intervals 
until  scavenging  is  terminated  at  the  point  where  the  sum  of  the  increments  of  mass  inflow  equal 
the  mass  capacity  of  the  combustor  at  scavenging  pressure  and  temperature,  that  is 


v 

2 AmiN  = 
x 


PsCAV  v0 
P Tscav 


(A4.16) 


The  combustor  temperature,  t0,  at  the  end  of  the  last  scavenging  time  interval  can  then  be 
compared  with  the  trial  value  of  rov  selected  at  the  start  of  the  blowdown  analysis.  If  there  is  a large 
discrepancy  between  these  two  values  of  temperature,  then  the  value  of  tov  can  be  modified 
and  the  cycle  calculations  repeated  until  the  agreement  is  acceptable. 


A4.3  Estimation  of  a Mean  Turbine  Efficiency 

The  mean  efficiency  over  the  whole  cycle  can  be  evaluated  from 

Y 

2 AW'act 

riTisis  = y (A4. 1 7) 

ZAIfis 

w 
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TABLE  I 


Details  of  Blade  Characteristics  at  Reference  Diameter  for  Single-Stage  Turbine 


Details  of  Stator  Row 

Inlet  Blade  Angle.  {3, 

0° 

Outlet  Blade  Angle.  (3,  (nominal) 

-60° 

Inlet  Gas  Angle,  a, 

0° 

Outlet  Gas  Angle,  a, 

-60° 

Blade  Pitch/Blade  Chord,  s/c 

0-75 

Maximum  Blade  Thickness/Blade  Chord.  /m»x/c 

0-200 

Radial  Tip  Clearance/Annulus  Height,  k/h 

0-010 

Annulus  Height/Blade  Chord,  hjc 

1-2 

Details  of  Rotor  Row 

Inlet  Blade  Angle.  |34 

0° 

Outlet  Blade  Angle.  (nominal) 

-60° 

Inlet  Gas  Angle.  a4 

Variable 

Outlet  Gas  Angle.  as 

-60° 

Blade  Pitch/Blade  Chord,  s/c 

0-75 

Maximum  Blade  Thickness/Blade  Chord.  tm»x/c 

0-200 

Radial  Tip  Clearance/Annulus  Height,  k/h 

0-020 

Annulus  Height/Blade  Chord,  h/c 

1-2 

TABLE  II 


Design  Point  Flow  Conditions  and  Flow  Areas  for  Single-Stage  Turbine 


Station 

Number 

B 

t 

K 

P 

kPa 

P 

kPa 

V 

m/s 

M 

a 

o 

P 

o 

(nom.) 

Aa 

mm2 

mm2 

A 

mm2 

I 

600  0 

589-9 

620-6 

584-9 

142-3 

0-292 

0-0 

0 

2607 

2607 

2607 

600-0 

559-7 

603-7 

473-3 

284-6 

0-600 

-60-0 

-60 

3057 

1529 

1529 

600-0 

559-7 

603-7 

473-3 

284-6 

0-600 

-60-0 

3057 

1529 

■ 

569-8 

559-7 

503-7 

473-3 

142-3 

0-300 

0-0 

0 

3057 

3057 

3057 

569-8 

529-5 

485-6 

375-7 

284-6 

0-617 

-60-0 

-60 

3644 

1822 

1822 

569-8 

529-5 

485-6 

375-7 

284-6 

0-617 

-60-0 

3644 

1822 

H 

539-6 

529-5 

401-3 

375-7 

142-3 

0-309 

0 0 

3644 

3644 

m = 1 • 283  kg/s 
U = 246-5  m/s 


TtisiO)  = 0- 1296 
Tt46(0>  = 0- 1647 


t.max  = 75-7%  (Determined  Graphically)  tjmax  = 75-41  % (Determined  Graphically) 

vorr  = 0-655  (Determined  Graphically)  vopt  = 0-614  (Determined  Graphically) 


TABLE  IV 

Details  of  Blade  Characteristics  at  Reference  Diameter  for  Two-Stage  Turbines 


Blade  Geometrical  Characteristic 

Turbine  Build  Number 

1 

2 

3 

4 

First  Stator 

Inlet  Blade  Angle.  p, 

0 

0 

0 

0 

Outlet  Blade  Angle.  (nominal) 

-60 

-50 

-70 

-60 

Inlet  Gas  Angle,  a, 

0 

0 

0 

0 

Outlet  Gas  Angle.  a2 

-60 

-50 

-70 

-60 

Blade  Pitch/Blade  Chord,  s/c 

0-75 

0-75 

0-75 

0-75 

Maximum  Blade  Thickness/Blade  Chord.  rmax/c 

0-200 

0-200 

0-200 

0-200 

Radial  Tip  Clearance/ Annulus  Height,  k/h 

0-010 

0-010 

0-010 

0-010 

Annulus  Height/Blade  Chord,  h/c 

1-2 

1*2 

1-2 

1-2 

First  Rotor 

Inlet  Blade  Angle.  jSt 

0 

0 

0 

30 

Outlet  Blade  Angle.  (nominal) 

-60 

-50 

-70 

-60 

Inlet  Gas  Angle.  a4  (design  value) 

o 

0 

o 

30 

Outlet  Gas  Angle.  a5 

-60 

-50 

-70 

-60 

Blade  Pitch/Blade  Chord,  s/c 

0-75 

0-75 

0-75 

0-75 

Maximum  Blade  Thickness/Blade  Chord.  /m»x/c 

0-200 

0-200 

0-200 

0-200 

Radial  Tip  Clearance/ Annulus  Height,  k/h 

0-020 

0-020 

0-020 

0-020 

Annulus  Height/Blade  Chord,  h/c 

1-2 

1-2 

1-2 

1-2 

Second  Stator 

Inlet  Blade  Angle.  J37 

0 

0 

0 

30 

Outlet  Blade  Angle.  [3,  (nominal) 

-60 

-50 

-70 

-60 

Inlet  Gas  Angle.  a7  (design  value) 

0 

0 

0 

30 

Outlet  Gas  Angle.  a8 

-60 

-50 

-70 

-60 

Blade  Pitch/Blade  Chord,  s/c 

0-75 

0-75 

0-75 

0-75 

Maximum  Blade  Thickness/Blade  Chord.  rm»x/c 

0-200 

0-200 

0-200 

0-200 

Radial  Tip  Clearance/Annulus  Height,  k/h 

0-010 

0-010 

0-010 

0-010 

Annulus  Height/Blade  Chord,  h/c 

1-2 

1-2 

1-2 

1-2 

Second  Rotor 

Inlet  Blade  Angle.  p,0 

0 

0 

0 

30 

Outlet  Blade  Angle.  pn  (nominal) 

-60 

-50 

-70 

-60 

Inlet  Gas  Angle.  a10  (design  value) 

0 

0 

0 

30 

Outlet  Gas  Angle.  au 

-60 

-50 

-70 

-60 

Blade  Pitch/Blade  Chord,  s/c 

0-75 

0-75 

0-75 

0-75 

Maximum  Blade  Thickness/Blade  Chord.  tm»x/e 

0-200 

0-200 

0-200 

0-200 

Radial  Tip  Clearance/ Annulus  Height,  k/h 

0-020 

0-020 

0-020 

0-020 

Annulus  Height/Blade  Chord,  h/c 

1-2 

1-2 

1-2 

1-2 

\ TABLE  V \ 

gn' Point  Flow  Conditions  and  Flow  Areas  for  the  Fow  Two-Stage  'Nirfoines 


tiH 


Turbine  Build  Number  1 


fi 

O 

(nom.) 


0 2607 

-60  3057 

3057 
0 3057 

-60  3644 

3644 
0 3644 

-60  4351 

4351 
0 4351 

-60  5296 


m = 1 • 283  kg/s 

TtisiO)  = 0- 1296 

U = 246  -5  m/s 

Tt4«(0i  = 0- 1647 

Tt7»(0>  = 0- 1296 

Ttioi2(<»  = 0- 1647 

Turbine  Build  Number  2 

600-0 

OO 

620  6 

546-9 

206-8 

0-429 

0-0 

0 

1882 

1882 

600-0 

548-5 

603-2 

440-7 

321-7 

0-686 

-50-0 

-50 

2214 

1423 

600-0 

548-5 

603-2 

440-7 

321-7 

0-686 

-50-0 

2214 

1423 

569-8 

548-5 

503-5 

440-7 

206-8 

0-441 

0-0 

0 

2214 

2214 

569-8 

518-3 

485-6 

348-5 

321-7 

0-706 

-50-0 

-50 

2645 

1700 

569-8 

518-3 

485-6 

348-5 

321-7 

0-706 

-50-0 

2645 

1700 

539-6 

518-3 

401-2 

348-5 

206-8 

0-453 

0-0 

0 

2645 

2645 

539-6 

488-1 

388-9 

273-8 

321-7 

0-727 

-50-0 

-50 

3171 

2038 

539-6 

488-1 

388-9 

273-8 

321-7 

0-727 

-50-0 

3171 

2038 

509-3 

488-1 

317-9 

273-8 

206-8 

0-467 

0-0 

0 

3171 

3171 

509-3 

457-8 

305-5 

210-4 

321-7 

0-751 

-50-0 

-50 

3871 

2489 

509-3 

457-8 

305-5 

210-4 

321-7 

0-751 

-50-0 

3871 

2489 

479-1 

457-8 

246-6 

210-4 

206-8 

0-482 

0-0 

3871 

3871 

m = 1-283  kg/s 
U = 246-5  m/s 


y-ns,o,  = 0- 1062 
JW,  =0-1303 
Tt7»(0i  = 0- 1062 

Ttioi2(Oi  = 0-1303 


TABLE  V — continued 

Design  Point  Flow  Conditions  and  Flow  Areas  for  the  Four  Two-Stage  Turbines 


Turbine  Build  Number  3 


Station 

Number 


p 

P 

kPa 

kPa 

620-6 

606-2 

601-6 

489-8 

601-6 

489-8 

502-1 

489-8 

481-1 

387-4 

481-1 

387-4 

397-6 

387-4 

384-2 

305-4 

384-2 

305-4 

314-0 

305-4 

P 

o 

(nom.) 


Ak  A n A 
mm2  mm2  mm 


241  -9  234-8  89-7 10-205 


0 4032 

-70  4737 


-70  5670 

5670 
0 5670 

-70  6785 


-70  8297 

8297 


m = 1 • 283  kg/s 

Tti3(0)  = 0- 1695 

U = 246  5 m/s 

Tt4««>>  = 0 2226 

1t79(0)  = 0- 1695 

Ftioi2(0i  = 0-2226 

Turbine  Build  Number  4 

\ 

600  0 I 588 
600  0 
600  0 
569-8 
569  8 
569-8 
539-6 
539-6 
539-6 
509-3 
509-3 
509-3 
479-1  |464 


71620-6 
601-8 
601-8 
502-2 
469  5 
469-5 
388-0 
366-5 
366-5 
299  6 
278-2 
278-2 
-0  1 224-6 


580-5 
457-1  : 
457-1 
457-1 
351-3 
351-3 
351-3 
269  6 : 
269-6 
269-6 
200-8 
200-8 
200-8 


0-0 
-60-0 
\60-0 
ifcl-O 
-60  0 
-60-0 
30  0 
-60-0 
-60  0 
30-0 
-60  0 
-60-0 
30-0 


0 2472 

-60  2958 

2958 
30  2958 

-60  3639 

. 3639 
30  3639 

-60  4469 

. 4469 
30  \4469 
-60  3633 

5633 


rh  = !•  283  kg/s 

Ttis,  o,  = 0-1296 

\ 

U = 174-3  m/s 

)W>  = 0-2753 

\ 

Ft7*(0)  = 0-2209 

\ 

Ftioi2(0)  = 0-2753 

\ 

* 

\ , 
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246-5m/s.  p13  = 101  -4  kPa.  Mean  Turbine  Efficiency:  ifnsi3  = 73-59%.  Mean  Turbine  Inlet  Temperature:  T,  (mass)  = 720-0  K.  ^(timk)  = 705- 0 K. 


Station  numbers  (numbers  in 
brackets  apply  to  second  stage) 


Sign  convention  for 
gas  and  blade  angles 


FIG.  1 DIAGRAMMATIC  REPRESENTATION  OF  TURBINE  STATOR  AND  ROTOR 
TOGETHER  WITH  VELOCITY  DIAGRAMS 


c Blade  chord 

s Blade  pitch 

t Maximum  blade  thickness 

max 

t Blade  trailing-edge  thickness 

i Incidence  angle  of  flow 

a Gas  flow  angle 

0 Blade  angle 


FIG.  2 TURBINE  BLADE  NOMENCLATURE 


Selected  blade  shape 


T6  profile  on  parabolic 
camber  line  - 1 /c  = 0.200 

mix 


Design-point  velocity  diagrams  as  selected  in  Section  2.1.1 
together  with  design-point  velocities  as  determined  in  Section  2.1 .3 


Design-point  velocity  diagrams 


i 


FIG.  3 SELECTED  BLADE  SHAPE  AND  DESIGN-POINT  VELOCITY  DIAGRAMS 
FOR  THE  SINGLE-STAGE  TURBINE 


Turbine  loss  characteristics  determined 
by  revised  Ainley  correlation 
Turbine  loss  characteristics  determined 
by  Modified  correlation 


OjO  0.1  02  03  0A  05  0.6  0.7  0.8  0.9  1.0  1.1 

Velocity  ratio  I>T1T13  fT1s13 


FIG.  7 EFFICIENCY-VELOCITY  RATIO  RELATIONSHIPS  - BUILDS  1 TO  4;  DESIGN  SPEED; 
LOSSES  DETERMINED  BY  BOTH  REVISED  AINLEY  AND  MODIFIED  CORRELATIONS 
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Percentage  of  cycle  time 

FIG.  9 COMBUSTOR  PRESSURE,  TEMPERATURE  AND  MASS  THROUGHOUT  PULSATING-FLOW  CYCLE 
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UNITED  KINGDOM 

Australian  Defence  Science  and  Technical  Representative  146 
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Short  Brothers  & Harland  183 
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Professor  G.  K.  Batchelor  187 

Sir  William  Hawthorne  188 
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